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ABSTRACT 

We study the efficiency of pair production in polar caps of young pulsars under a variety of conditions 
to estimate the maximum possible multiplicity of pair plasma in pulsar magnetospheres. We develop 
a semi-analytic model for calculation of cascade multiplicity which allows efficient exploration of the 
parameter space and corroborate it with direct numerical simulations. Pair creation processes are 
considered separately from particle acceleration in order to assess different factors affecting cascade 
efficiency, with acceleration of primary particles described by recent self-consistent non-stationary 
model of pair cascades. We argue that the most efficient cascades operate in the curvature ra¬ 
diation/synchrotron regime, the maximum multiplicity of pair plasma in pulsar magnetospheres is 
^ few X 10®. The multiplicity of pair plasma in magnetospheres of young energetic pulsars weakly 
depends on the strength of the magnetic field and the radius of curvature of magnetic field lines and 
has a stronger dependence on pulsar inclination angle. This result questions assumptions about very 
high pair plasma multiplicity in theories of pulsar wind nebulae. 

Keywords: acceleration of particles — plasmas — pulsars: general — stars: neutron 


1. INTRODUCTION 

The idea that production of electron-positron pairs in 
magnetospheres of rotation-powered pulsars is intimately 
connected with their activity had been proposed by Stur- 
rock (1971) only a few years after the discovery of pulsars 
(Hewish et al. 1968). Since then it has become an inte¬ 
gral part of the standard pulsar model and today, there 
is little doubt that an active rotationally powered pulsar 
produces electron-positron plasma. Although the pulsar 
emission mechanism(s) is still not yet identified, there is 
strong empirical evidence that pulsars stop emitting ra¬ 
dio waves when pair formation ceases - the threshold for 
pair formation roughly corresponds to the “death line” 
in pulsar parameter space, what was already noted by 
Sturrock (1971). Furthermore, the narrow peaks in many 
pulsar high-energy light curves (Abdo et al. 2010) require 
pervasive screening of the whole magnetosphere by pair 
plasma, except in narrow accelerator gaps (e.g. Watters 
et al. 2009; Pierbattista et al. 2015). Understanding pair 
plasma generation in pulsar magnetospheres is there¬ 
fore of crucial importance for developing pulsar emission 
models. 

In the standard pulsar model, initially proposed by 
Goldreich & Julian (1969) and Sturrock (1971), the mag¬ 
netosphere is filled with dense pair plasma which screens 
the accelerating electric field everywhere except some 
small zones which are responsible for particle accelera¬ 
tion and emission. Pair plasma is primarily created via 
conversion of 7 rays in the strong magnetic field near the 
polar caps (PCs). Pair production in the PCs is a “cor¬ 
nerstone” of the standard model - without dense plasma 
produced at the PCs, at the base of open magnetic field 
lines, the magnetosphere would have large volumes with 
unscreened electric field, as pair creation in e.g. outer 
gaps (Cheng et al. 1976) cannot screen the electric field 
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over the rest of the magnetosphere. 

Charge starvation (Arons & Scharlemann 1979) or vac¬ 
uum gaps (Ruderman & Sutherland 1975) at the polar 
cap, when the number density of charged particles is not 
enough to screen the electric field, leads to formation 
of accelerating zone(s). Some charged particles enter 
this zone, are accelerated to very high energies and emit 
gamma-rays which are absorbed in the ultrastrong mag¬ 
netic field, creating electron-positron pairs. The pairs, 
being relativistic, can also emit pair producing photons 
and so the avalanche develops until photons emitted by 
the last generation of pairs can no longer produce pairs 
and escape the magnetosphere. 

The pair plasma created by pulsars flows out of the 
magnetosphere along open magnetic field lines and pro¬ 
vides the radiating particles for the surrounding Pulsar 
Wind Nebulae (PWNe). Models of PWNe depend (at 
least) on the density of the plasma, what produces the 
observed synchrotron and inverse Compton emission. Es¬ 
timates of the pair multiplicity (the number of pairs pro¬ 
duced by each primary accelerated particle) needed to 
account for the emission from the Crab pulsar wind neb¬ 
ula (PWN) range from about 10® — 10® (de Jager et al. 
1996) up to lO’^ (e.g. Bucciantini et al. 2011); for the Vela 
PWN the multiplicity is estimated to be about 10® (de 
Jager 2007). PWNe therefore give the most compelling 
evidence for pair production and pair cascades in at least 
young energetic pulsar magnetospheres. 

Although PWNe are observed only around young pul¬ 
sars (< a few times 10'^ yrs), evidence for pairs, at least 
for high plasma densities larger than those provided by 
primary particles, can also be found in older pulsars. 
Synchrotron absorption models for the eclipse in the dou¬ 
ble pulsar system PSR J0737-3039 (Arons et al. 2005; 
Lyutikov 2004) require a pair multiplicity of around 10® 
for the recycled 22 ms pulsar in that system. 

The cascade process in pulsar polar caps has been 
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the subject of extensive numerical as well as analytical 
studies (e.g. Daugherty & Harding 1982; Gurevich & Is¬ 
tomin 1985; Zhang & Harding 2000; Hibschman & Arons 
2001a, b; Medin & Lai 2010). The pair plasma multi¬ 
plicity obtained in these studies was significantly lower 
than estimates of pair plasma multiplicity in PWNe, as 
it did not exceed ~fewxl0^. Most of those works con¬ 
sidered pair creation together with the particle accel¬ 
eration which makes these analyses dependent on the 
acceleration model considered. These studies also as¬ 
sumed steady, time-independent acceleration of the pri¬ 
mary particles. However, recent studies by Timokhin 
(2010); Timokhin & Arons (2013) have found that pulsar 
polar cap pair cascades are not time-steady in the gen¬ 
eral case of arbitrary current, particularly those required 
by global magnetosphere models (e.g Contopoulos et al. 
1999; Timokhin 2006; Spitkovsky 2006; Kalapotharakos 
& Contopoulos 2009). 

In this paper we study the question of what is the max¬ 
imum pair multiplicity achievable in pulsar polar cap cas¬ 
cades and under which circumstances is it achieved. In 
contrast to previous pair cascade studies, we take a mul¬ 
tistep approach. We consider the physical processes in 
pair cascades and particle acceleration models separately 
in order to clearly set apart different factors influencing 
the efficiency of pair cascades. We first assess how each 
of the microscopic processes affects the final multiplic¬ 
ity and the pulsar parameter ranges that result in the 
largest possible pair multiplicity. Then, we employ the 
most recent model of non steady-state particle accelera¬ 
tion in pulsar polar caps and derive a simple analytical 
estimate for the maximum energy of particles accelerated 
in a non-stationary cascade. One on the most important 
results of our study is a strong upper limit on pair plasma 
multiplicity in pulsars. 

We limit ourselves to the case of cascades at the polar 
caps of young pulsars as from previous theoretical stud¬ 
ies of polar cap cascades, such pulsars are expected to 
be the most efficient pair producers. We rely on results 
of previous cascade studies in our choice of the specific 
cascade process, namely cascades initiated by curvature 
radiation of primary particles. 

The plan of the paper is as follows. In § 2 we briefly 
discuss the efficiency of cascades in general and give an 
overview of the most efficient cascade process in pulsar 
polar caps. In subsequent sections we consider in detail 
all physical processes in such cascades. §§3-7, the largest 
part of this paper, are devoted to development of a simple 
semi-analytical model for estimation of pair production 
efficiency in polar caps of young pulsars. This model al¬ 
lows efficient exploration of pulsar parameter space and 
helps to clarify the main factors affecting cascade multi¬ 
plicity. Then, in §8 we use direct numerical simulations 
of PC cascades in a Crab-like pulsar to show that pre¬ 
dictions of our semi-analytical model are indeed correct. 
We discuss uncertainties of current pulsar models in § 10 
and summarize and discuss our findings and their impli¬ 
cations in the Discussion. 

2 . PHYSICS OF POLAR CAP CASCADES: AN 
OVERVIEW 

An electron-positron cascade can be thought of as a 
process of splitting the energy of primary particles into 
the energies of secondary particles. The maximum mul¬ 


tiplicity, the number of secondary particles for each pri¬ 
mary particle, of an ideal cascade initiated by a single 
primary particle with energy ep would be 

Kmax ^ 2 . (1) 

^ 7 , esc 

where esc is the maximum energy of photons escaping 
from the cascade (or the minimum energy of pair pro¬ 
ducing photons). Not all of the primary particle energy 
goes into pair producing photons, and pairs created in 
the cascade do not radiate all of their energy into the 
next generation of pair producing photons. Hence, Kmax 
is only the theoretical upper limit on the multiplicity of 
a real cascade. The total pair yield of the cascade is a 
combination of four factors: (a) number of primary par¬ 
ticles, (b) initial energy of primary particle - the higher 
the energy the more pairs can be produced, (c) threshold 
for pair formation - the lower the threshold the higher 
the multiplicity, (d) efficiency of splitting the energy of 
primary particles into pairs - the higher the fraction of 
particle energy going into pair creating photons, as op¬ 
posed to the final kinetic energy of the particles and pho¬ 
tons below the pair formation threshold, the higher the 
multiplicity. 

In young, fast rotating pulsars, the electric field in the 
polar cap acceleration zones is strong and primary par¬ 
ticles can be accelerated up to very high Lorentz factors, 
7 > 10^. At these energies the most effective radiation 
process is curvature radiation (CR). CR efficiency grows 
rapidly with the particle energy and for young pulsars 
becomes the dominant emission mechanism for primary 
particles. For secondary particles, which are substan¬ 
tially less energetic than the primary ones, the primary 
way to create pair producing photons is via synchrotron 
radiation. In § 9 we argue that although another possible 
emission mechanism for secondary particles - Resonant 
Inverse Compton Scattering (RICS) of soft X-ray pho¬ 
tons emitted by the NS - can generate pair producing 
photons in pulsars with magnetic field > 3 x 10^^ G, 
that channel never becomes the dominant one for pair 
production in young pulsars, at best resulting in pair 
multiplicity comparable to the one of the synchrotron 
channel. Hence, in young pulsars which are expected to 
have the highest multiplicity pair plasmas, the polar cap 
cascades should operate primarily in the CR-synchrotron 
regime - all known studies of cascades in pulsar polar cap 
agree on this point. In this paper we study in detail CR- 
synchrotron cascades. The resulting multiplicities will 
be good estimates for a wide range of pulsar parame¬ 
ters, however, as we consider only synchrotron radiation 
of secondary particles, for pulsars with magnetic field 
> 3 X 10^^ G our analysis might underestimate the cas¬ 
cade multiplicity by a factor of ~ 2, see § 9. 

Fig. 1 gives a schematic overview of electron-positron 
cascade development in polar cap regions of young pul¬ 
sars. Shown are the first two generations in a cascade 
initiated by a primary electron. Primary electrons emit 
CR photons (almost) tangent to the magnetic held lines; 
primary electrons and CR photons are generation 0 par¬ 
ticles in our notation. Magnetic held lines are curved 
and the angle between the photon momentum and the 
magnetic held grows as the photon propagates further 
from the emission point. When this angle becomes large 
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Figure 1. Schematic representation of electron-positron cascade 
in the polar cap of a young pulsar, see text for description. 

enough, photons are absorbed and each photon creates 
an electron-positron pair - generation 1 electron and 
positron. The pair momentum is directed along the mo¬ 
mentum of the parent photon and at the moment of cre¬ 
ation, the particles have non-zero momentum perpendic¬ 
ular to the magnetic field. They radiate this perpendic¬ 
ular momentum almost instantaneously via synchrotron 
radiation and then move along magnetic field lines. Al¬ 
though these secondary particles are relativistic, their 
energy is much lower than that of the primary electron 
and their curvature photons cannot create pairs. After 
the emission of synchrotron photons, secondary particles 
(generation 1 and higher) no longer contribute to cas¬ 
cade development. Generation 1 photons (synchrotron 
photons produced by the generation 1 particles) are also 
emitted (almost) tangent to the magnetic field line - as 
the secondary particles are relativistic - and propagate 
some distance before acquiring the necessary angle to 
the magnetic field and creating generation 2 pairs. These 
pairs in their turn radiate their perpendicular momentum 
via synchrotron radiation, emitting generation 2 photons. 
The cascade initiated by a single CR photon stops at a 
generation where the energy of synchrotron photons falls 
below e^^esc- 

Only primary particles emit pair producing photons 
as they move along the field lines; all secondary par¬ 
ticles emit pair producing photons just after their cre¬ 


ation. The cascade development can be thus divided 
into two parts: (i) primary particles emit CR photons 
as they move along magnetic field lines and (ii) each CR 
photon gives rise to a synchrotron cascade, when syn¬ 
chrotron photons create a successive generation of pairs 
which emit the next generation synchrotron photons at 
the moment of creation. This division goes between gen¬ 
eration 0 and all subsequent cascade generations. 

In the following sections we analyze all four factors 
regulating the yield of CR-synchrotron cascades listed at 
the beginning of this section (in reverse order, from d to 
a). In §§3-7 we develop a simple semi-analytical mod¬ 
els for calculation of multiplicity of strong polar caps 
cascades initiated by a primary electron which we then 
compare with detailed numerical computation described 
in §8. We start with gamma-ray absorption in a strong 
magnetic field in §3; then we discuss the efficiency of the 
synchrotron cascade in §4. Curvature radiation is consid¬ 
ered in §5 and in §5.3 we discuss the multiplicity of a CR- 
synchrotron cascade initiated by a single particle of given 
energy. This covers items c and d from the list of fac¬ 
tors affecting cascades efficiency. In §6 we address item b 
from the list - energy of primary particles initiating the 
cascade. §§7 and 8 are devoted to the total multiplicity 
of a CR-synchrotron cascade when particle acceleration 
is taken into account - in §7 we present results for cascade 
multiplicities from a semi-analytical model and in §8 we 
present results of detailed numerical simulations of cas¬ 
cades in a Crab-like pulsar. In §9 we discuss the role of 
RICS in polar caps cascades and argue that considering 
only CR-synchrotron cascades gives us an adequate esti¬ 
mate of the pair multiplicity in young pulsars. Finally, 
in § 10 , we address item a - the mean flux of primary par¬ 
ticles and the total yield of a CR-synchrotron cascade in 
an energetic pulsar. We argue that this is the most im¬ 
portant factor regulating pair yield in energetic pulsars. 
Despite the uncertainty in determining the mean flux of 
primary particles, we can set a rather strict upper limit 
on pair multiplicity in pulsars. 

3. PHOTON ABSORPTION IN THE MAGNETIC 
FIELD 

3.1. Opacity for 7 — H pair production 

The opacity for single photon pair production in strong 
magnetic field is (Erber 1966) 

aB(e 7 , V’) = 0.23 ^ 5 sini/'exp (2) 

where b = B/Bg is the local magnetic field strength 
B normalized to the critical quantum magnetic field 
Bq = ejafX^ = 4.41 x 10^^ G, is the angle between 
the photon momentum and the local magnetic field, 
a/ = !he « 1/137 is the fine structure constant, and 
Ac = h/mc = 3.86 x 10“^^ cm is the reduced Compton 
wavelength. The parameter y is defined as 

y=ie^ 6 sinV', (3) 

where is the photon energy in units of irieC^. For con¬ 
venience from here on, all particle and photon energies 
will be quoted in terms of rUeC^. The optical depths for 
pair creation by a high energy photon in a strong mag- 
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Pc=10^cm 
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Figure 2. Contour plot of 1/xa as a function of the logarithms of the magnetic field strength B in Gauss, and photon energy normalized 
to the electron rest energy, for three values of the radius of curvature of magnetic field lines pc = 10®, 10^, 10®cm. Ixa values shown on 
this plot are calculated from eq. (7) and are not corrected for the kinematic threshold (see text). 


Pc=10®cm 


Pc=10^cm 


Pc=10^cm 



Figure 3. Contour plot of 1/xa corrected for kinematic threshold according to eq. (9). These values of 1/xa are used in all calculations 
within semi-analytic model. Notations are the same as in Fig. 2. 


netic field after propagating distance I is 

r(e^,/)= / aB{e.y,'ip{x))dx, (4) 

^0 

where integration is along the photon’s trajectory. 

Expression (2) is accurate if the magnetic field is small 
compared to the critical field Bg, b < 0.2 suffices, and 
if sin ip > 2 so that the created pair is in a high Lan¬ 
dau level; pair production into low Landau levels and 
for higher magnetic fields has been discussed in Daugh¬ 
erty & Harding (1983) and Harding et al. (1997). For 
most pulsars, the magnetic field in polar cap regions is 
smaller that Bg/i. In this paper we study cascades in 
pulsars with “normal” magnetic fields and so we neglect 
high-held effects. 

Expression (2) becomes inaccurate when pairs are cre¬ 
ated at low Landau levels, near the pair formation thresh¬ 
old, it overestimates the opacity and even formally allows 


pair formation below the kinematic threshold sin ip = 
2. However, pairs created at low Landau levels will not 
give rise to strong cascades, as their perpendicular en¬ 
ergy will be too low to emit pair-producing synchrotron 
photons. Hence, for the case of strong cascades, accurate 
treatment of pair formation near the kinematic threshold 
is not necessary. Throughout this paper we use Erber’s 
approximation (2) in all our analytical calculations but 
introduce a cut off at the kinematic threshold e-y sin ip = 2 
as described at the end of this section, thus taking into 
account the cessation of pair formation below the the 
kinematic threshold. 

If the mean free path (mfp) of a pair-producing pho¬ 
ton ly is comparable to or larger than the characteristic 
scale of the magnetic field variation Lb, this photon will 
not initiate a strong cascade with the same emission pro¬ 
cesses by which it was produced. The reason for this is 
as follows. The opacity for 'yB pair production exponen¬ 
tially depends on the magnetic held strength and photon 
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energy via x- The energy of the next generation photon 
will be smaller than that of the primary one, and, be¬ 
cause the primary photon has already traveled the dis¬ 
tance over which the magnetic field has substantially de¬ 
clined, the magnetic field along the next generation pho¬ 
ton’s trajectory will be substantially weaker than that 
along the primary photon’s trajectory^. The next gener¬ 
ation photon’s mfp will be much larger than than that of 
the primary photon, and, even if this secondary photon 
will be absorbed, it will be the last cascade generation. 
Hence, in a strong cascade, for all but the last gener¬ 
ation photons, Ij Lb- A reasonable estimate for Lb 
would be the distance of the order of the NS radius R^s as 
any global NS magnetic field decays with the distance as 
6 > 3. Very localized, sun spot like magnetic 
fields, are in our opinion of no importance for the gen¬ 
eral pulsar case as the probability of such a “spot” to lie 
at the polar cap should be rather low, i.e. most pulsars 
should be able to produce plasma in a more or less regular 
magnetic field. A dipole field, d = 3, is often considered 
as a reasonable assumption for a general pulsar model. 
Pure dipole field, however, seems to be a too idealized 
approximation, as even if the NS field is a pure dipole, it 
will be slightly disturbed by the currents flowing in the 
magnetosphere. In general, near dipole magnetic fields 
with different curvatures of magnetic field lines should 
be examined in cascade models. 

We consider strong cascades with large multiplicities, 
where, as argued above, photons propagate distances 
much shorter than the characteristic scale of the mag¬ 
netic field variation, so we assume that in the region 
where most of the pairs are produced the magnetic field 
is constant. The radius of curvature of magnetic field 
lines Pc is not smaller than Lb, and as I-y ^ Lb the 
angle ip is always small, the approximation sin^ « i/; is 
very accurate. For photons emitted tangent to the mag¬ 
netic field line, dx = Pcdip. In our approximation both b 
and Pc are constants. From eq. (3) we have ip = 2x/cyb, 
and substituting it into eq. (4) we can express the optical 
depth T to pair production as an integral over x 


more efficient numerical codes than direct integration of 
eq. (5). 

As the optical depth to pair formation grows exponen¬ 
tially with X (and distance), for analytical estimates it is 
reasonable to assume that all photons are absorbed when 
they have traveled the distance ^ such that T{l,y) = 1 . 
We denote the value of x when the optical depth reaches 
I as Xa : 

Xa-T(Xa) = l- (7) 

In all our computations we will use Xa as the value of the 
parameter x at the point of the photon’s absorption. Xa 
is a solution of the non-linear equation (7) with t given 
by eq. ( 6 ). Because of the exponential dependence of r 
on 1 /x it is to be expected that 1 /xa should have a close 
to linear dependence on logarithms of e-y, b, and pc- We 
solved equation (7) numerically for different values of ey, 
b, and pc and, indeed, the inverse quantity I /xa depends 
almost linearly on log Cy, log b, and log pc in a wide range 
of these parameters. Making a modest size table of 1/xo 
values one can later use piecewise interpolation to find a 
particular value of Xa quite accurately. 

In Fig. 2 we plot contours of 1/xa as functions of 
log(e.y) and log(H) for three different values for the ra¬ 
dius of curvature of the magnetic field lines. We want to 
point out that 1 /xa differs from the value 1 /xa = 15 used 
by Ruderman & Sutherland (1975), especially for higher 
energy photons. Although this difference is only a factor 
of a few, as we will point out later, the cascade efficiency 
in each generation depends on the corresponding value of 
Xa, and in a strong cascade with several generations, the 
estimate for the final multiplicity will be substantially 
affected by the value of Xa- 

Erber’s expression is not applicable for pairs produced 
at low Landau levels as it overestimates the opacity and, 
formally, solutions for Xa obtained from eq. (7) allows 
pair creation even for e.ySini/) <2. In our analyti¬ 
cal treatment we introduce a limit on Xa to correct for 
the kinematic threshold. For photons above kinematic 
threshold from eq. (3) it follows that 


xiCy, /) 




e^yb 


Xexp 


w) ''*■ 


( 5 ) 


where A^- = 0.92ay/Ao « 1.74 x 10® cm“^. 

Integrating eq. (5) over x by parts two times we can 
get an expression for r in terms of elementary functions 
and the exponential integral function Ei: 


t{x) = A^ 




1 4 \ 



( 6 ) 


Ei(z) is a widely used special function, defined as Ei(z) = 
— J^^exp{—t)/t dt. There are efficient numerical algo¬ 
rithms for its calculation implemented in many numer¬ 
ical libraries and scientific software tools; using eq. ( 6 ) 
for calculation of the optical depths will result in much 


^ In the polar cap, photons are emitted by ultrarelativistic parti¬ 
cles moving along magnetic field lines. At emission points photons 
are almost tangent to field lines and the differences in initial photon 
pitch angles for different generations can be neglected. 


Xa > b , (8) 

In all our analytical calculation we get Xa from eq. (7) 
and then use the maximum value of Xa and b : 

Xa = max(xa, 5) ■ (9) 

In Fig. 3 we plot contours of 1/xa which incorporate cor¬ 
rections to Xa due to the kinematic threshold according 
to eq. (9). It is clear from the plot that this correction af¬ 
fects only cases with high magnetic field and low particle 
energies. 


3.2. Energy of escaping photons 

As discussed above, photons escaping the cascade are 
those whose mfp ly is larger than the distance of signifi¬ 
cant magnetic field attenuation Lb- The formal criteria 
we use for calculating the energy of escaping photons 
esc ^ 7 (^- 7 ,esc) '^esc-^NSi ^esc dimeusionlcss pa¬ 
rameter quantifying the escaping distance in units of i?Ns. 
The photon mfp is ly = Pcipa and expressing the angle 
between the photon momentum and the magnetic field 
at the point of absorption ipa through Xa we get a (non- 
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log Gy,esc 


created pairs. 


13.0 



6.0 6.5 7.0 7.5 8.0 

log Pc 


Figure 4. Energy of escaping photons: contours of log esc ^is a 
function of logarithms of the radius of curvature of magnetic field 
lines pc in cm and magnetic field strength B in Gauss for Sesc = 1- 


linear) equation for e.^,esc 


^^ 7 , esc 


= 2 


Pc Xa(^ 7 , esc; Pc) 

iSC.RNS b 


( 10 ) 


The non-linearity in this equation is due to dependence 
of Xa on e.y,esc; b, and Pc- Using an interpolation table 
for 1/xa this equation is very easy to solve numerically 
for all reasonable values of physical parameters involved. 

In Fig. 4 we plot energy of escaping photons, log esc 
as a function of the radius of curvature of magnetic field 
lines Pc and magnetic field strength B for Sesc = 1- For 
smaller values of Sesc the whole plot would move to the 
left. This figure shows (an obvious) trend that for higher 
magnetic field and smaller radii of curvature, the energy 
of escaping photons is lower, which allows for more cas¬ 
cade generations and larger multiplicity. The break in 
contour lines around logB ~ 12.4 is due to the kine¬ 
matic threshold. 

For a given geometry of the magnetic field a more ac¬ 
curate analytical expression for the energy of escaping 
photons can be derived using an exact expression for p 
along photon’s travel path by expanding the integral in 
eq. (4) around its upper endpoint, as was done by Hib- 
schman & Arons (2001b). However here we try to explore 
different possible magnetic field conhguration - explor¬ 
ing parameter space in pc - and our simple estimate is 
accurate to a factor of a few. 


4. SYNCHROTRON CASCADE 

In this section we discuss the synchrotron cascade, 
where most of the electron-positron pairs are created. 
The synchrotron cascade is the part of the whole cas¬ 
cade that is initiated by generation 0 photons. In the 
synchrotron cascade each generation’s primary photon is 
divided into many (lower energy) next generation’s pair- 
producing photons by synchrotron radiation of freshly 


4.1. Fraction of parent photon energy remaining 
in the cascade 

A high energy photon, when absorbed in the magnetic 
field, produces an electron and a positron; the total en¬ 
ergy of these particles is equal to the energy of the pho¬ 
ton. Particle momenta just after production have pitch 
angles equal to pa, when pair production takes place well 
above threshold. When pairs are created at high Landau 
levels, as is the case in strong cascades, relativistic par¬ 
ticles have non-zero pitch angles and they radiate their 
perpendicular energy via synchrotron radiation; in su¬ 
perstrong magnetic fields, this happens almost instanta¬ 
neously. The component of particle momentum parallel 
to the magnetic field is unaffected by synchrotron radi¬ 
ation and so the final Lorentz factor of the particle p 
will be 

e±,F = (l~/ 3 y) ~ e±, I [l + (lAo e±, i)^] ^ (H) 


where /3|| = uy/c is particle velocity along the magnetic 
field line and e±^i is the initial Lorentz factor of the par¬ 
ticle right after creation. If the photon absorption hap¬ 
pens at Xa < 1 “ which is indeed the case near pulsar 
polar caps, see §3.1 - with a high degree of accuracy we 
can assume that the energy of the photon is equally di¬ 
vided between the electron and positron (e.g. Daugherty 
& Harding 1983). Expressing pa in eq. (11) through Xa, 
pa = 2xap'yb, and the initial particle energy through the 
photon energy e.^, e±,i = e^/2 we get e±,F as a function 
of Xa and b 


e±.F = 



( 12 ) 


The fraction of photon energy radiated as syn¬ 
chrotron photons, which is going into subsequent pair 
creation, is 


Csyn — 


2(e±.i - e±.F) 


= 1 - 



(13) 


Because of the kinematic threshold (8) the minimum 
value of this fraction is Csynlxa=b — 0.292, i.e. formally^ 
at least ~ 30% of absorbed photon energy will go into 
synchrotron radiation of created pairs. In Fig. 5 we plot 
the fraction of pair-producing photon energy radiated as 
synchrotron photons by freshly created pairs given by 
eq. (13). Contours of Csyn are plotted as functions of the 
photon energy e-y and magnetic field strength H, the ra¬ 
dius of curvature was assumed to be pc = lO^cm. The 
dependence of Csyn on pc (via Xa) is very weak, and Fig. 5 
is a good representation of how Csyn depends on and 
B for any pc of interest. 

It is evident from Fig. 5 that for higher magnetic field 
strengths, B > 3 x lO^^G, and lower energies of parent 
photons, a progressively smaller fraction of the parent 
photon energy goes into synchrotron photons; the rest 
remains in the kinetic energy of the created pairs mov¬ 
ing along magnetic field lines. The portion of the parent 


^ as we mentioned above, the physics of near-threshold pair for¬ 
mation is more complicated and our simplified treatment is less 
accurate in this regime. 
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lOgCy 

Figure 5. Fraction of the parent photon energy radiated as syn¬ 
chrotron photons by freshly created pairs: contours of Csyn a 
function of logarithms of the parent photon energy and mag¬ 
netic field strength B in Gauss for pc = lO^cm. 


lOgCy 

Figure 6. Number of synchrotron photons with the characteristic 
energy e^^syn emitted at each 7 —> e+e“ conversion: contours of 
Usyn as a function of logarithms of the parent photon energy 
and magnetic field strength B in Gauss for pc = lO^cm. 


photon energy energy left in kinetic energy of pairs does 
not go into production of next generation pairs but is 
“lost” from the synchrotron cascade^. The reason for 
this is that for higher magnetic held strengths pairs are 
created when the photon has a smaller pitch angle ijja, 
so that a smaller fraction of the photon energy goes into 
perpendicular pair energy, and hence, a smaller fraction 
of the photon energy is emitted and remains in the cas¬ 
cade. 


4.2. Number of secondary photons 

At each pair creation event the parent photon is ef¬ 
fectively transformed into an electron-positron pair and 
lower energy synchrotron photons. Those photons be¬ 
come parent photons for the next cascade generation or 
escape the magnetosphere, terminating the cascade. 

The characteristic energy of synchrotron photons emit¬ 
ted by a newly created particle in terms of quantities used 
in this paper is given by 

3 

£ 7,8711 = (14) 

The number of synchrotron photons with the character¬ 
istic energy e-y^ syn ~ these photons carry most of the en¬ 
ergy of synchrotron radiation - emitted at each event 
of conversion of a parent photon with energy Cj into an 
electron-positron pair is 


..., ^7 

'^syn — 

£7,3711 


4 Csyn 
3 Xa 


(15) 


In eq. (15) both ^syn and Xa are functions of B, e-y, and 
Pc- In Fig. 6 we plot contours of risyn as functions of the 


^ the kinetic energy of pairs might be tapped by RIGS cascade 
branches, see § 9 


photon energy and magnetic field strength, the radius 
of curvature was assumed to be pc = lO^cm. Two clear 
trends are visible on this plot: the lower the energy of the 
primary photon the larger the number of secondary syn¬ 
chrotron photons produced at each conversion event, and 
the higher the magnetic held the smaller is the number 
of synchrotron photons. The hrst one is a general trend 
of emission processes when higher energy particles emit 
less photons which, however, have larger energies. The 
second trend is due to the suppression of the synchrotron 
cascade discussed above, in §4.1. 


4.3. Multiplicity of synchrotron cascade 

The generation i -|- 1 cascade photon is a synchrotron 
photon emitted at the event of pair creation by a photon 

(i) 

of generation i. Expressing b, if a and through Xa and 
from eq. (14) we get for the characteristic energy of 
the next generation photon 


^7 


3 

4 


Xa£ 


(i) 

7 ’ 


(16) 


where and are energies of i’th and (i + l)’th 

generation photons. The photon energy degrades with 
each successive generation of the cascade. This degrada¬ 
tion accelerates as the cascade proceeds through genera¬ 
tions because with the decrease of the photon energy Xa 
increases, see Fig. 2. 

The photon mfp in a constant magnetic held goes as 
ly (X Xal^-y- For photon energies < lO'^, when Xa ^ 10, 
the mfp of successive generations increases by at least an 
order of magnitude in each successive generation; there¬ 
fore the crudeness of our approximation for estimating 
the energy of escaping photons, eq. ( 10 ), could affect 
only the last cascade generation. 

The rapid energy degradation results in a rather small 
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number of generations in polar cap cascades as the energy 
of pair-producing photons rapidly reaches the threshold 
energy. On the other hand, the number of emitted syn¬ 
chrotron photons increases with the decrease of the en¬ 
ergy of the parent photon, see Fig. 6 ; consequently syn¬ 
chrotron cascades can have quite large multiplicities de¬ 
spite a small number of generations. 

Each cascade photon creates 2 particles at the mo¬ 
ment of absorption, and particles are produced until 
£ 7 ^^^ > e..^ esc- The total number of particles gener¬ 
ated in synchrotron cascades initiated by a primary pho¬ 
ton can be calculated by summation over all generations 
of the number of pairs produced in each cascade gen¬ 
eration. The algorithm for this calculation is shown in 
Appendix A, algorithm 1 . We will use this algorithm 
for calculation of the polar cap cascade multiplicity after 
we discuss curvature radiation, the radiation mechanism 
responsible for generating the primary photons for syn¬ 
chrotron cascades in young energetic pulsars. 

Finally we wish to point out how the magnetic field 
strength affects the multiplicity of the synchrotron cas¬ 
cade. From discussions presented in §§3.2, 4.1 it is ev¬ 
ident that, for the same energy of the primary photon, 
a higher magnetic field strength results in (i) reducing 
the final energy of escaping photons, therefore increasing 
cascade efhciency, and at the same time (ii) in suppres¬ 
sion of cascade efficiency by forcing a smaller fraction of 
the photon energy to go into perpendicular energy of the 
created pairs and so short-circuiting the cascade. There¬ 
fore, there should be a “sweet spot” in magnetic field 
strength where the synchrotron cascade is most efficient. 


5. CURVATURE RADIATION 

In this section we discuss how primary particles emit 
photons which “launch” the synchrotron cascade. As we 
discussed above the most efficient process for supplying 
the primary (generation 0 ) photons in young energetic 
pulsars is curvature radiation. 


5.1. Fraction of the primary particle energy going into 
the cascade 


^CR 


8.0 


7.5 


O 

Q. 
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o 
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loge® 

Figure 7. Fraction of the primary particle energy emitted as CR 
photons over the distance scr = 1 - contours of ^cr as a function 
of logarithms of the initial particle energy and the radius of 
curvature of magnetic field lines pc in cm. 



Solving eq. (18) we get for the particle energy after it 
travels the distance s from the injection point 


e±(.s) = 4 


(el)3 


-1/3 


(19) 


(see also Harding 1981), where s is normalized to i?NS: 
is the initial particle energy; constant H is defined as 
= (2/3)i?Nsi'e ~ 1-88 X 10“^cm^, where r^ = e^ImgC^ 
is the classical electron radius. 

The fraction ^cr of the initial particle energy lost to 
curvature radiation after a particle has traveled distance 
Scr is 

Ccr(scr) = 1-'#^^. (20) 

4 


Ultrarelativistic particles moving along curved mag¬ 
netic field lines emit electromagnetic radiation with 
power (Jackson 1975) 


Pr.n = 


3 nrigC^ \p, 


4 - 


(17) 


where -fcR is normalized to /sec, e± is the particle’s 
energy normalized to m(,c^; we do not distinguish be¬ 
tween electrons and positrons. Since we are treating the 
pair generation problem separately from the problem of 
particle acceleration, we consider cascades produced by 
particles injected in a region with screened electric field, 
so that particles are not accelerated and only lose en¬ 
ergy to CR. The particle’s energy decreases with time 
according to the equation of motion 


de± 

dt 


= -Pr 


(18) 


If the energy of these CR photons goes into creation of 
electron-positron pairs, Ccr gives the efficiency of the CR 
part of the full cascade. The electric field in the ac¬ 
celeration zone transforms electromagnetic energy into 
particle’s kinetic energy, which is then radiated as pair 
producing photons. Only the photon’s energy can be di¬ 
vided in chunks carried by a large number of pairs. The 
cascade will have high efficiency if (i) primary particles 
have high energy, (ii) emit most of their energy as pho¬ 
tons, and (iii) inject these photons in the region where 
the synchrotron cascade can work effectively, i.e. in a 
region close to the NS which is smaller than the charac¬ 
teristic scale of magnetic field variation L^- 
In Fig. 7 we plot the fraction of the primary particle 
energy emitted as CR photons after the particle has trav¬ 
eled distance Scr = 1. Shown are contours of Ccr(scr) as 
a function of the initial particle energy and the radius 
of curvature of magnetic field lines pc- For smaller val¬ 
ues of Scr the whole plot would move to the right. CR is 
most efficient in transferring particle energy into the cas- 
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log e±,th 



log Pc 

Figure 8. Critical particle energy above which it can emit pair- 
producing photons via curvature radiation: contours of loge^ 
as a function of logarithms of the radius of curvature of magnetic 
field lines pc in cm and magnetic field strength B in Gauss for 
-5esc = 1- 

cade in the parameter space corresponding to the lower 
right triangular region of Fig. 7. Going from the upper 
left (smaller e^, larger pc) to the lower right (larger 
smaller pc) on this plot, not only the particle energy in¬ 
creases but also the fraction of the energy which can go 
into the cascade. 

For a certain range of and pc the value of the energy 
put into the cascade by the primary particle grows faster 
than the energy of that particle. For more or less regular 
global magnetic fields, with pc > lO^cm, the transition 
between effective and ineffective CR cascades occurs at 
particle energies e± ~ 10^, and the efficiency is very sen¬ 
sitive to the particle energy. In this parameter range 
even a modest increase of the primary particle energy 
can result in a large boost of the cascade multiplicity. 

5.2. Energy of CR photons and critical particle energy 

The characteristic energy of CR photons emitted by a 
particle with the energy e± is (Jackson 1975) 

e7.cR= ^ —e±"«5.8x 103p-i (21) 

^ Pc 

where pc ,7 = Pc/lO^cm and e ±_7 = e±/10^. The number 
of CR photons emitted by the particle while traveling 
distance ds normalized to i?Ns is 

dn-cR ^ Rns Pcr (22) 

ds c OR 

Each CR photon above the pair formation threshold will 
be a primary photon for the synchrotron cascade dis¬ 
cussed in §4. The critical energy at which primary par¬ 
ticles can produce pair-creating CR photons can be cal¬ 
culated by equating e-y^cR given by eq. (21) to the escape 
photon energy e-^^esc from eq. (10). In Fig. 8 we plot 


the critical particle energy which could initiate pair pro¬ 
duction with CR photons e±^th for Sesc = 1- Shown are 
contours of loge±_th as a function of the radius of curva¬ 
ture of magnetic field lines pc and magnetic field strength 
B. Primary particles should have energies e± > 10® to 
be able to initiate pair production via CR. 


5.3. Multiplicity of CR-synchrotron cascade 

The multiplicity of the the CR-synchrotron cascade - 
the total number of particles produced by a single pri¬ 
mary electron or positron accelerated in the gap - can 
be computed by multiplying the number of CR photons 
ricR, eq. (22), by the number of particles produced in 
the synchrotron cascade initiated by these photons risyn, 
eq. (15), and integrating it over the distance where CR 
can initiate a cascade. 


^CR-syn 



d ricR 


ds . 


(23) 


The actual algorithm we use to compute the total mul¬ 
tiplicity is the Algorithm 2 from Appendix A. In Fig. 9 
we plot AtcR-syn ss a function of the initial particle en¬ 
ergy loge± and magnetic field strength logR for three 
different radii of curvature of magnetic field lines pc = 
10®, 10^, 10® cm assuming Scr = Sesc = 1- 

It is evident from these plots that in a dipolar magnetic 
field, with pc « 10®cm, the maximum achievable multi¬ 
plicity is AvoR-syn ~ few X 10"*^ even in cascades initiated 
by extremely energetic primary particles. If the radius 
of curvature is an order of magnitude less, a rather high 
multiplicity Kcnsyn ^ 10® could be achieved in polar cap 
cascades for magnetic field strength B ^ 10^^ G and 
particle energies e± > lO’’, parameters quite realistic for 
young pulsars. For strongly non-dipolar magnetic field, 
with Pc « 10®cm the multiplicity can be another order 
of magnitude higher ACcR-syn few x 10®. 

The properties of CR do not depend on the strength 
of the magnetic field, therefore the effect of the magnetic 
field strength on the multiplicity of CR-synchrotron cas¬ 
cade is due to the synchrotron cascade and how many 
CR photons pair produce. As discussed at the end of 
§ 4.3 there should be an optimum magnetic field strength 
where the multiplicity is the highest; the multiplicity de¬ 
creases both for higher (due to larger energy left in par¬ 
ticle motion along magnetic field lines) as well as lower 
magnetic field (due to increase of the energy of escap¬ 
ing photons). This trend is clearly visible on all plots 
of Fig. 9. The highest cascade multiplicity is achieved 
for magnetic field around B ~ 10^^ G. The value of this 
optimum magnetic field grows slightly with increase of 
Pc but it stays around ^ few x 10^^ G even for dipolar 
magnetic fields. This is noteworthy in view of the fact 
that B ~ 10^^ G is the typical value of magnetic field 
strength for normal pulsars. For any given energy of the 
primary particle the decrease of the cascade multiplicity 
towards stronger magnetic fields is faster than for weaker 
fields. 

The dependence of cascade multiplicity on the initial 
energy of the primary particle is non-linear. Let us con¬ 
sider what happens when the initial particle energy 
goes from the highest to the lowest value (horizontal di¬ 
rection in plots of Fig. 9). For the highest values of 
multiplicity decreases uniformly, but then it drops by an 
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Figure 9. Multiplicity of CR-synchrotron cascade: contours of log/^cR-syn as a function of logarithms of the primary particle energy 
and magnetic field strength B in Gauss for three values of the radius of curvature of magnetic field lines pc = 10®, 10^, lO^cm. Assumed 
values for characteristic lengths: scr = 5esc = 1- 
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Figure 10. Efficiency of CR-synchrotron cascade as given by eq. (24): contours of CcR-syn as a function of logarithms of the primary particle 
energy and magnetic field strength B in Gauss for three values of the radius of curvature of magnetic field lines pc = 10®, 10^, 10®cm. 
Assumed values for characteristic lengths: scr = Sesc = 1- 



order of magnitude in a rather small range of after 
about half a decade of values the multiplicity drops 
again to 1, when no particles can be produced. The first 
effect is due to the decrease of the efficiency of CR, dis¬ 
cussed in §5.1. For lower initial energies of the primary 
particle there is less energy available to create pairs, not 
only because the particle energy is smaller, but also due 
to smaller efficiency of CR in producing photons - the 
primary particle keeps most of its energy, depositing only 
a small fraction of it in the cascade zone. The drop in the 
efficiency of CR Cor) where it becomes less than 10% (see 
Fig. 7) manifests in a rapid decrease of cascade multiplic¬ 
ity - by an order of magnitude - on all plots of Fig. 9 for 
magnetic field strengths where the maximum multiplic¬ 
ity is achieved. This drop in KoR-syn is less pronounced 
for both higher and lower magnetic field strengths due 
to lower efficiency of the cascade discussed above. The 


second drop in AtcR-syni towards 1, is due to the threshold 
in pair formation - for those particle energies CR pho¬ 
tons have too low an energy to initiate a cascade. The 
blue region on the plots of Fig. 9 show the parameter 
space where no particles can be produced by the CR- 
synchrotron cascade. It does not mean, however, that 
no pairs can be produced in the polar cap cascades for 
such primary particle energies. Instead of CR, the pri¬ 
mary photons for synchrotron cascade will be produced 
by inverse Compton scattering of thermal photons emit¬ 
ted by the NS, however, those primary photons will have 
much lower energies and multiplicities of such cascades 
will be quite low (see Harding & Muslimov 2002). 

We find it also instructive to compare the multiplicity 
of the CR-synchrotron cascade with the theoretical upper 
limit on cascade multiplicity Kmax, given by eq. (1) in §2. 




















Pair multiplicity in young pulsars 


11 


The ratio 

(24) 

^max 

can be considered as the efficiency of splitting the en¬ 
ergy of the primary particle into pairs. In Fig. 10 we 
plot CcR-syn for the same values of parameters as KcR-syn 
in Fig. 9. Despite lower multiplicity for smaller values 
of B, the cascade efficiency is higher, i.e. more of the 
initial energy of the primary particle goes into pair for¬ 
mation as opposed to the energy of escaping photons and 
kinetic energy of pairs and primary particles. This trend 
is discussed in §4.1, 5.1. It is interesting to note that for 
B < few X 10^^ G the cascade can be quite efficient in 
splitting a noticeable fraction of primary particle energy 
into pairs. For magnetic field ~ 10^^ G the fraction of 
the primary particle’s energy going into pair production 
saturates at ~ 10%; this is the limiting efficiency of the 
highest possible multiplicity cascade in a typical pulsar. 
The dependence of CcR-syn on is similar to the depen¬ 
dence of KcR-syn On e. 

6. PARTIGLE AGGELERATION 
6.1. Overview of particle acceleration regimes 

In this section we will get an estimate for the energy 
of primary cascade particles. In the following discus¬ 
sion we will rely on results of self-consistent modeling of 
pair cascades by Timokhin (2010) [TIO] and Timokhin 
& Arons (2013) [TA13]. First, we give a brief overview 
of how particle acceleration proceeds according to these 
simulations. 

Whether and how the pair formation along given mag¬ 
netic field lines occurs depends on the ratio jm/jcj of the 
current density required to support the twist of magnetic 
field lines in the pulsar magnetosphere (e.g. Timokhin 
2006; Bai & Spitkovsky 2010), jm = (c/47r) V x B, to the 
local GJ current density, jcj = 77gjC, where tjqj = B/Pc 
is the GJ charge density. 

For the Ruderman & Sutherland (1975) cascade model, 
where particles cannot be extracted from the NS sur¬ 
face, effective particle acceleration and pair formation 
is possible for almost all values of jm/joj (TIO). In the 
space charge limited flow regime, hrst discussed by Arons 
& Scharlemann (1979), pair formation is not possible if 
0 < Jm/jcj < 1, but is possible for all other values of 
Jm/joj (TA13). Pair formation is always non-stationary: 
an active phase when particles are accelerated to ultra- 
relativistic energies and give rise to electron-positron cas¬ 
cades - with a burst of pair formation - is followed by 
a quiet phase when recently generated dense electron- 
positron plasma screens the electric field everywhere. 

As pair plasma leaves the active region, it flows into 
the magnetosphere, and later into the pulsar wind. When 
the density of the pair plasma drops below the minimum 
density necessary for screening of the electric field a gap 
appears - a charge starved region where the electric field 
is very strong, of the order of the vacuum electric field. 
To screen the electric field, the plasma density must be 
high enough to provide both the GJ charge density tjqj 
and the imposed current density The transition be¬ 
tween the region(s) still filled with plasma (hereafter we 
call it “plasma tail”) and the gap is sharp - plasma still 
capable of screening the electric held moves in bulk, its 
density is close to the critical density. The plasma den¬ 


sity drops abruptly at the gap boundary; within the gap 
the particle number density is much smaller than the 
critical density (cf. Fig. 22, 23 in § 10). The motion 
of the boundary between the plasma Hlled region and 
the gap sets the gap’s growth rate. Some particles enter 
the gap and are accelerated to ultrarelativistic energies. 
The gap grows until the energy of particles accelerated 
there becomes sufficient to start electron-positron cas¬ 
cades and the cycle repeats. The gap is not station¬ 
ary; in almost all cases it moves as a whole - after its 
growth has been terminated by newly created pairs its 
upper and low boundaries are moving in the same di¬ 
rection. The gap can move, keeping its size for a long 
time, or it might disappear rather quickly. The details 
of gap behavior - where the gap appears, what direction 
it moves, how fast it disappears - depend on the ratio 
jm/jcj- However, despite these differences, the way in 
which the highest energy particles are accelerated is very 
similar in any regime which allows pair formation stud¬ 
ied in TIO, TA13. Namely, the size of the charge starved 
region grows as the tail of pair plasma moves, the bulk 
velocity of the tail v sets the rate of the gap expansion. 
Particles entering the gap from the tail are accelerated 
in a larger and larger gap, until they are able to produce 
pair-producing photons. The place where these photons 
are absorbed and produce pairs is the other boundary of 
the gap. 

6.2. Energy of primary particles 

In this section we obtain a quantitative estimate for 
the maximum energy of accelerated particles using as 
an example the case of the Ruderman-Sutherland (RS) 
cascade, when particles cannot be supplied from the sur¬ 
face of the NS. As we mentioned before, gap formation 
and particle acceleration in the space charge limited flow 
regime, when pair creation is allowed, is very similar to 
the RS case and estimates for particle energies obtained 
in this section are applicable for the space charge limited 
flow with jm/jcj > 1 and jm/jci < 0. The GJ charge 
density is positive and we consider the case when the ra¬ 
tio jm/jcj > 0. In Fig. 11 we show a schematic picture of 
how particles are accelerated in this cascade. On the top 
of each figure we show the electric field in the accelerat¬ 
ing region and on the bottom a schematic representation 
of plasma motion in and around the gap; plot (a) corre¬ 
sponds to the time when the electric field screening has 
just started, plot (b) shows a well-developed gap moving 
into the magnetosphere. These schematic plots illustrate 
results of actual simulations of RS cascades shown in 
Fig. 3,4,11 and 13 in TIO. 

At the beginning of the burst of pair formation, the 
gap appears at the NS surface and its upper boundary is 
the “tail” of plasma left from the previous burst of pair 
formation, where the particle number density is still high 
enough to screen the electric field (region I in Fig. 11). 
Electrons and positrons in this tail are trapped in elec¬ 
trostatic oscillations and the bulk velocity of this tail v 
is sub-relativistic, but for large current densities (around 
or greater than Jgj) it is quite close to c. Electrons from 
this tail which get to the gap boundary are pulled into 
the gap and are accelerated toward the NS. As the tail 
moves, the gap grows; the current and charge density in 
the gap is due to the flux of electrons from the tail and 
so it remains constant within the gap. The gap growth 
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Figure 11. Schematic representation of gap formation and evolution for cascades in Ruderman-Sutherland regime with jm/jcj > 0. See 
text for explanation. 


is stopped when electrons reach an energy high enough 
to produce pair-creating photons. This first-generation 
of pairs start screening the electric field - electrons move 
toward the NS and positrons are accelerated toward the 
magnetosphere and start producing pair creating pho¬ 
tons as well (region II in Fig. 11(a)). Because those 
positrons are ultrarelativistic they practically co-move 
with the photons and so new pairs are injected close to 
their parent particles making a blob of pair plasma mov¬ 
ing into the magnetosphere (region Ila in Fig. 11(b)). 
This blob is the lower boundary of the accelerating gap, 
and the gap exists until this blob catches up with the tail 
from the previous pair formation cycle. For large current 
densities this can take a while as v is close to c. Plasma 
leaking from the blob forms the new tail (region lib on 
Fig. 11(b)). 

The evolution of the electric field in any given point x 
and moment of time t is given by (see e.g. eq. 1 in TA13) 

QE 

— (a;,t) = -47r(j(a:,t) - jm) = -47rJ, (25) 

j is the actual current density along a given magnetic 
field line and jm is the current density imposed by the 
magnetosphere. The difference j = j — in the gap 
remains constant. When the upper boundary of the gap 
moves with the constant speed v this equation can be 
integrated to get the electric field in the gap 

E{x, t) = E{xq, to) + 47rJ ——- + Attj (t - to) (26) 

V 


introduce defined as 




Jgj 



(28) 


where is the GJ current density in an aligned rotator 


Jg,i = ’IgjC = 


B 


(29) 


is a factor which shows how stronger/weaker the elec¬ 
tric field in the gap is compared to the situation of a static 
vacuum gap in an (anti-)aligned rotator, like the one con¬ 
sidered by Ruderman & Sutherland (1975). v ~ const is 
a good approximation to the numerical results and so 
~ const. In cascades along magnetic held lines where 
jm is close to the local value of Jgj in an aligned rotator 
~ 2, for the same situation in a pulsar with inclination 
angle of 60°, Ci ~ 1- 

For energy losses dominated by curvature radiation, 
free acceleration is a good approximation for B > 10^^ G 
(see Appendix B). If radiation losses are negligible, the 
particle’s equation of motion is 


dp 

dt 


= —eE , 


(30) 


where p = mcj is particle’s momentum. In terms of the 
distance traveled by the particle in the gap l± = c(t — to) 
with E given by eq. (27), particle equation of motion can 
be written as 


dp 



(31) 


Where E{xo, to) is the electric held within the gap at the 
moment to at the point Xo- If we assume that the gap 
boundary is at Xo at the moment to, then E{xo,to) = 0. 

Electrons enter the gap from above and are quickly 
accelerated by the strong electric held and move with 
relativistic speed practically from the moment they leave 
the plasma tail. If a particle enters the gap at to (in point 
Xo) its coordinate is x± = xo — c{t — to), substituting x± 
into eq. (26) the electric held seen by that particle is 
given by 


Ayr / C \ n 

E{xp) = —J (1 -h - ) (xo - x±) = 47r?7°j l± (27) 
c V V/ 


Integrating this equation and expressing jcj through pul¬ 
sar parameters we get for particle energy 




27r 

EqXcC 



(32) 


The distance the primary particle travels in the region 
of unscreened electric held Igap - the size of the gap as 
seen by the moving particle - is the sum of the distance 
the particle travels before emitting pair-producing pho¬ 
tons terminating the gap ^±,gap and the distance these 
protons travel until the absorption point gap 

^gap = ^±,gap + ^7, gap ■ 


In the last step in eq. (27) we denote with l± = \x± — For any given particle, the larger the distance Zj. the par- 

XqI the distance traveled by the particle in the gap and tide travels to the emission point, the higher the particle 
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energy and the energy of CR photons it emits, and so 
the smaller is the distance traveled by the photon until 
the absorption point l-y. The distance the particle trav¬ 
els in the gap Igap is the minimum value oi I = + Ij 

because once the first pairs are injected the avalanche of 
pair creation will lead to screening of the electric field. 
The photon mean free path l-y can be estimated from 
eq. (3) as 

lj=2xa^. (34) 

Photon energy depends on the particle energy e± 
which depends on according to eq. (32), and so 
is a function of Zj.. Zgap can be found by minimizing 
I = using as an independent variable: Zgap is 

the value of Z which satisfies dl/dl% = 0, gap and Zj. gap 
are the values of l-y and Z^ where Z reaches its minimal 
value. Using eq. (33) we can write an equation for gap 


dl 

~dl 


7. gap 


= - 1 . 


±,gap 


(35) 



If the photon energy depends on Zj_ as e-y oc (Zj_)“, then 
eq. (35) is reduced to 

^±,gap ~ ^ (36) 


Figure 12. Primary particle energy: contours of loge^^acc as a 
function of logarithms of the radius of curvature of magnetic field 
lines pc in cm and magnetic field strength B in Gauss. We used 
the following values for gap parameters P = 33 ms, ^-j = 2 and 
Xa = 1/7. 


where gap is expressed in terms of Z^_ using eq. (34). 
Zgap is then given by 

I — ('37') 

‘gap” ‘i.gap- 


particles accelerated in the gap"^ 


acc — 


A/T „4I7 p-l!7 Ta-117 

\ Xlc) 


49 f irBa 
18 


5 X 10^ xT if ■ (41) 


The final energy of the primary particle is given by 
eq. (32) with l± = Zgap. Please note that because the 
gap moves, the actual size of the gap (see Fig. 11(a)) is 

^gap “ 4" ^gap ■ (38) 

The energy of the CR photons depends on the particle 
energy as (eq. 21), the particle energy depends on l± 
as Z^ (eq. 32), hence, €j oc (Z^)® and a = 6. Substituting 
expression for e± (eq. 32) into the expression for CR pho¬ 
ton energy Cj (eq. 21), the latter into eq. (34), and the 
resulting expression for into eq. (36) with a = 6 after 
algebraic transformations we get the following expression 
fO'' ^±,gap 


Z 


e 

± 


>gap 


1/7 

■ (39) 


The size of the gap as seen by the moving particle ac¬ 
cording to eq. (37) is 

Zgap O. 2 X 104 xTC^'p^J} cm , (40) 

where B 12 = R/IO^^G and pc .7 = Pc/lO^cm. Substitut¬ 
ing l± = Zgap into eq. (32) we get for the final energy of 


The dependence of the primary particles energy on pul¬ 
sar period P, filling factor and the strength of the 
magnetic field B is very weak, the only substantial de¬ 
pendence is on the radius of curvature of magnetic field 
lines. The reason for this is the strong dependence of 
CR photon energy on the energy of emitting particles 
Cj oc e^, eq. (21). Changes in the threshold energy of 
pair producing photons which can stop the gap growth 
cause only modest variation of the energy of primary 
particles, which explicitly depends on P, and B, but 
not on Pc; eq. (32). The energy of pair producing pho¬ 
tons sets the energy of accelerated particles, and in pul¬ 
sars with a strong accelerating electric field the gap will 
be smaller than in pulsars with a weaker accelerating 
field. In Fig. 12 we plot the energy of particles e±,acc 
accelerated in the gap of a pulsar with P = 33 ms as 
a function of the radius of curvature of magnetic field 
lines Pc and magnetic field strength B, assuming = 2 
and Xo = 1/7- The value Xa ~ 1/7 corresponds to 
Xo of CR photons emitted by relativistic particles with 

= 2.5 X 10^ in a magnetic field B = lO^^ G with 

Our expression for the energy of primary particles has the 
same dependence on pc, P and B as the expression for the poten¬ 
tial drop in the gap derived by Ruderman Sz Sutherland (1975), 
their eq. (23). This is to be expected as in both cases particles 
are accelerated by the electric field which grows linearly with the 
distance and the size of the gap is regulated by absorption on cur¬ 
vature photons in magnetic field. The difference is in the presence 
of factor and a different numerical factor. 
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Figure 13. PP diagram with the yellow area showing the 
range of parameters where approximation for particle acceler¬ 
ation used in this paper is applicable, see text for descrip¬ 
tion. Pulsars from ATNF catalog (Manchester et al. (2005), 
http://www.atnf.csiro.au/research/pulsar/psrcat) are shown by 
black dots, 7 -ray pulsars from the second Fermi catalog (Abdo 
et al. 2013) by red dots. 

Pc = 10^ cm, this is a good estimate for Xa in eq. (41) 
for young pulsars as the dependence on Xa is very weak. 
This plot clearly illustrates the dependence of e±_ acc on 
B and pc - weaker magnetic field and/or larger radius of 
curvature requires larger photon energies for terminating 
gap growth, and so the energy of the primary particles 
is larger. 

We derived eq. (41) under the following assumptions: 
(i) particles are accelerated freely, i.e. radiation reac¬ 
tion can be neglected, (ii) the length of the gap is much 
smaller that the polar cap radius, so that a one dimen¬ 
sional approximation can be used, (iii) the magnetic field 
is B < 0.2Bq « 8.8 x 10^^ G so that the opacity to 'jB 
pair creation is described by eq. (2). Constraints on the 
pulsar parameters (i) and (ii) are derived in Appendix B 
and Appendix C correspondingly. Plotted on the PP 
diagram, Fig. 13, these restrictions select the range of 
pulsar period and period derivatives - shown as yellow 
region - where all three assumptions are valid. In this 
figure, the one-dimensional approximation (ii) is valid to 
the left of the solid line, given by eq. (C5), the approxi¬ 
mation (i) of free acceleration above the dot-dashed line, 
given by eq. (B2), and pulsars with B < 0.2Bq are below 
the dotted line. We see that most of young normal pul¬ 
sars, including gamma-ray pulsars from the Fermi second 
pulsar catalog, fall in this range. 

7. CASCADE MULTIPLICITY PER PRIMARY 
PARTICLE: SEMI-ANALYTICAL MODEL 

We now combine the results of §5.3 concerning the cas¬ 
cade multiplicity for a fixed energy of the primary parti¬ 
cle, and of §6.2 concerning the energy of primary parti¬ 
cles accelerated in the gap with given parameters P, 
strength H, and pc- The multiplicity of CR-synchrotron 


cascades depends on the energy of the primary particle 
e±, magnetic field B, and radius of curvature of magnetic 
field lines pc- The only significant dependence of the en¬ 
ergy of accelerated particles e± in young pulsars is on the 
radius of curvature of magnetic field lines Pc] the depen¬ 
dence on P, and B is very weak, see eq. (41). There¬ 
fore, when particle acceleration is taken into account, the 
overall cascade multiplicity k can depend substantially 
only on B and pc- 

In CR-synchrotron cascade changes in e± and pc 
change k in opposite directions - for higher e± multi¬ 
plicity is higher, for larger pc multiplicity is lower. But 
the energy of the primary particles accelerated in polar 
caps of young pulsars e± is higher for larger radii of cur¬ 
vature of magnetic field lines pc- Hence, increasing pc 
lowers the cascade multiplicity for a fixed e±, but at the 
same time increases e±, which partially compensates the 
decrease of cascade multiplicity. Therefore, the final cas¬ 
cade multiplicity should have a rather weak dependence 
on Pc: which leaves the magnetic field strength B the 
only parameter significantly affecting the multiplicity of 
strong cascades in polar caps of young pulsars. 

In Fig. 14 we show the final multiplicity as a func¬ 
tion of magnetic field strength and radius of curvature 
of magnetic field lines. We present plots for three sets 
of parameters P and which differ by ^ an order of 
magnitude. As expected, both pulsar period as well as 
the filling parameter have a very small effect on the 
final multiplicity. For the range of pc and B in Fig. 14, 
the CR-synchrotron cascade has the highest efficiency - 
primary particles lose most of their energy as CR pho¬ 
tons within the distance Sor = Rns where the cascade is 
possible - and higher multiplicities cannot be achieved. 
It is easy to see that the cascade multiplicity is not very 
sensitive to pc and for typical pulsar magnetic field of 
IQi^ G is around 10^ 

It is a remarkable fact that the multiplicity of the most 
efficient cascades is sensitive only to the strength of the 
magnetic field. For fixed B the total pair yield - the 
total number of particles injected into the magnetosphere 
- depends then only on the flux of primary particles. 
This is true for pulsars where the accelerating potential is 
regulated by pair production and where cascade operate 
in CR-synchrotron regime have high efficiency. 

8. CASCADE MULTIPLICITY PER PRIMARY 
PARTICLE: NUMERICAL SIMULATIONS 

In §§3-7 we developed a semi-analytical model of 
strong polar cap cascades. In order to verify the key 
assumptions and conclusions of this model we have per¬ 
formed numerical simulations of time-dependent polar 
cap pair cascades. Because such numerical simulations 
are quite time consuming we limited ourselves to the case 
of a young pulsar with parameters similar to the Crab 
pulsar (P = 33 ms) and magnetic field strength B and 
radius of curvature of magnetic field lines pc having val¬ 
ues resulting in high multiplicity. Our goal was to show 
that the main assumptions and conclusions of our anal¬ 
ysis of polar cap cascades are realistic. More extensive 
self-consistent numerical studies of the polar cap cascades 
will be done in subsequent papers. 

As outlined in §6.1, particles are quickly accelerated in 
the gap which is much smaller than the typical distance 
over which the full cascade develops. The primary pair 
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Figure 14. Multiplicity of polar cap cascades: contours of log k as a function of logarithms of curvature of magnetic field lines pc in cm 
and magnetic field strength B in Gauss for three sets of the gap parameters (P[ms],^j): (33,2), (33,0.25), (330,2). In all cases we used 
Xa = 1/7 for calculation of the energy of primary particles. 


producing particles are moving most of the time in the 
region with screened electric field. If the primary parti¬ 
cle energies are known, the full cascade can be modeled 
using traditional Monte-Carlo techniques (Daugherty & 
Harding 1982); to obtain the energies of primary particles 
initiating the CR-synchrotron cascade a self-consistent 
model of the cascade (Timokhin 2010; Timokhin & Arons 
2013) is necessary. 

We have performed numerical simulations of time- 
dependent polar cap pair cascades in a two-step process. 
In the first step, we use a hybrid Particle-in-Cell/Monte 
Carlo (PIC/MC) code PAMINA (PIC And Monte-Carlo 
code for cascades IN Astrophysics) to simulate the initial 
self-consistent electric field generation, particle accelera¬ 
tion and electric field screening near the NS to obtain 
the distribution functions of the electrons and positrons 
in the acceleration/screening region. This code includes 
only CR of the particles and first generation of pairs 
needed to screen the gap, and so does not follow the 
full synchrotron cascade. The details of this code are de¬ 
scribed in Timokhin (2010); Timokhin & Arons (2013). 

In the second step, we use another code to simulate the 
full pair cascade in the pulsar dipole field above the PC, 
including both CR of primary particles and synchrotron 
radiation of pairs. This code, based on the calculation de¬ 
scribed in detail in Harding & Muslimov (2011) [HMll], 
is a Monte-Carlo simulation of the electron-positron pair 
cascade generated above a PC by accelerated particles 
in the region of screened electric field. Although HMll 
included a steady particle acceleration component, this 
component is not used in the present calculation. We 
therefore assume that the particles and the further pairs 
they create do not undergo any acceleration. This code 
takes as input the distribution functions of accelerated 
particles output by the time-dependent PAMINA code 
that are moving away from the NS surface and simu¬ 
lates the combined cascade from all of the particles. Al¬ 
though our setup is capable of calculating full cascades 
generated by primary particles with arbitrary distribu¬ 
tion functions, in the simulations described in this section 
we used a monochromatic injection of primary particles. 


On the one hand, the energy distribution of the most en¬ 
ergetic primary particles which produce the bulk of the 
pairs in many cases is close to monochromatic (see e.g. 
§10, Fig. 23), and on the other hand this enables us to 
compare numerical simulations with predictions of our 
semi-analytical theory. The energy of the primaries was 
calculated from the self-consistent model however. 

The MC code first follows the primary particle in dis¬ 
crete steps along the magnetic field line at magnetic co¬ 
latitude 9, starting from the location Xq and particle en¬ 
ergy at time tpeak at the peak of the pair production 
cycle in PAMINA code, computing its curvature radia¬ 
tion. The steps Ax are set to the minimum of a fraction 
0.1 of a NS radius and the distance over which the par¬ 
ticle would lose 1% of its energy to curvature radiation. 
Dividing the CR spectrum at each step into logarith¬ 
mic energy intervals, a representative photon from each 
energy interval is followed through the curved magnetic 
field until its point of pair production (determined as a 
random fraction of the mean-free path). The number of 
CR photons in each energy bin, noR, is determined by 
the energy loss rate and average energy in that bin. The 
pairs produced by the photon, or the escaping photon 
number, is then weighted by Rcr- The created pair is 
assumed to have the same direction and half the energy 
of the parent photon. Although the CR photons are ra¬ 
diated parallel to the magnetic field, they must acquire 
a finite angle to the field before producing a pair, so the 
created pairs have finite pitch angles at birth. Each mem¬ 
ber of the pair emits a sequence of cyclotron and/or syn¬ 
chrotron photons, starting from its initial Landau state 
until it reaches the ground state, assuming the position 
of the particle remains fixed (given the very rapid radi¬ 
ation rate). As described in HMll, when the pair Lan¬ 
dau state is larger than 20, the asymptotic form of the 
quantum synchrotron rate (Sokolov & Ternov 1968) is 
used to determine the photon emission energy and final 
Landau state. When the Landau state is below 20, the 
full QED cyclotron transition rate (Harding & Preece 
1987) is used. At large distances above the NS surface, 
when the magnetic field drops below 0.002Bq, we short- 
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cut the individual emission sequence and use an expres¬ 
sion for the spectrum of synchrotron emission for an elec¬ 
tron that loses all of its perpendicular energy (Tademaru 
1973). Each emitted photon is then propagated through 
the magnetic held from its emission point until it pair 
produces or escapes. The next generation of pairs are 
then followed through their synchrotron/cyclotron emis¬ 
sion sequence. By use of a recursive routine that is called 
upon the emission of each photon, we can follow an arbi¬ 
trary number of pair generations. The cascade continues 
until all photons from each branch have escaped. As each 
member of each created pair completes its synchrotron 
emission, its ground-state energy, position and genera¬ 
tion number are stored in a pair table. As each photon 
either pair produces or escapes, its energy, generation 
and position of pair creation or escape are stored in ta¬ 
bles for absorbed and escaping photons. The photons 
and pairs from all accelerated particles are summed to¬ 
gether to produce the complete cascade portrait at that 
time step. 

The NS magnetic field in the MC code described above 
is a distorted dipole with an azimuthal (</)) component 
which is off-set from the center of the NS. The magnetic 
field is given by 


B = Bo 



X 


er cos[0(l -I- a)] + ^ sin[0(l + a)] 

—e^^e{9 + sm6cos6)sm{(j)—(j)o) , 


(42) 


where Bq is the surface magnetic field strength at the 
magnetic pole, r is the radial coordinate, a = e cos{(p — 
(j)o) is the parameter characterizing the distortion of po¬ 
lar held lines, and (j)o is the magnetic azimuthal angle 
dehning the meridional plane of the offset PC. The pa¬ 
rameter s sets the magnitude and the parameter (j)Q sets 
the direction of asymmetry of this azimuthal component. 
Setting e = 0 gives a pure dipole held structure, while a 
non-zero value of e produces an effective offset of the PC 
from the dipole axis in the direction specihed by (pQ. For 
non-zero e, the radius of curvature of the magnetic held 
lines is smaller than dipole in the direction of the offset 
and larger than dipole in the direction opposite to the 
direction of offset. This particular parametric form for 
the magnetic held was used in simulations of stationary 
cascades in HMll, it was chosen to account for distortion 
of the shape of the PC caused by currents howing in pul¬ 
sar magnetosphere. Such azimuthal asymmetries in the 
near-surface magnetic held is caused by the sweepback of 
the held lines near the light cylinder due to retardation 
(e.g. Dyks & Harding 2004) and currents (e.g. Timokhin 
2006; Bai & Spitkovsky 2010; Kalapotharakos et al. 2014) 
or, additionally, by asymmetric currents in the NS. 

We did not perform a systematic study of all parame¬ 
ter space with our numerical simulations, which will be 
done elsewhere, with our numerical simulations we test 
assumptions and predictions of our semi-analytic cascade 
model. Any form of magnetic held with adjustable radius 
of curvature of magnetic held lines would serve our pur¬ 
poses, but using the magnetic held given by eq. (42) al¬ 
lows comparison with the most recent simulations in the 


frame of the previous-generation cascade models HMll. 
We explored cases of pair cascades both for pure dipole 
helds and for azimuthally distorted helds. Multiplicities 
obtained from the numerical simulations agree reason¬ 
ably well with the semi-analytic model, within a factor 
of a few. As an example we describe in detail results of 
particular simulations with pulsar parameters yielding 
high multiplicity for a cascade at the peak of the pair 
creation cycle. The magnetic held is B = 10^^ G and is 
moderately distorted, with the offset e = 0.4 resulting in 
the radius of curvature of magnetic held lines near the 
NS Pc = 8.8 X 10® cm. The initial energy of primary 
particles from PAMINA simulations is = 2.3 x 10^. 

In Figures 15 and 16 we show a set of hve plots which 
we call “cascade portraits”. These plots illustrate dif¬ 
ferent aspects of the cascade development by showing 
moments of the photon or particle distribution function 
f{x,E). The top panel shows the number of particles 
produced in each energy and distance bin 

dN{x, E) = f{x, E) dx dE, (43) 


as a 2D color map. The number of particles is color 
coded in logarithmic scale according to the color bar on 
the right. The middle left plot shows the energy distri¬ 
bution E f{xi, E) of particles produced at four distances 
Xi, * = 1... 4; different line styles correspond to different 
distances according to the left plot legend. These spec¬ 
tra are essentially cross-sections of the map of dN{x, E) 
(multiplied by particle energy) along four lines shown in 
the plot for dN{x, E). The bottom left plot shows the en¬ 
ergy distribution of all particles produced in the cascade 
{dN(x, E) integrated along x direction and multiplied by 

E) 

Ef{E) = E dxf{x,E). (44) 

Jo 

In this and the following plots, colored lines show contri¬ 
butions of different cascade generations; lines are color- 
coded according to the right plot legend. The middle 
right plot shows the differential pair production rate - 
number of particles produced in a distance bin {dN(x, E) 
integrated along the E direction) 

dN{x)=dx / f{x,E)dE. (45) 

Jo 

The bottom right plot shows the cumulative pair produc¬ 
tion rate - total number of particles produced up to the 
distance x 


pExna.x 

N{<x)= dx f{x,E)dE. (46) 

Jo Jo 

Figure 15 shows the cascade portrait of the pairs (we do 
not differentiate between electrons and positrons). The 
cascade extends to about ^ 6i?NS) where it dies out com¬ 
pletely as the magnetic held strength and the acceler¬ 
ated particle energy decrease. The pair number grows 
very quickly, within a few tenths of a NS radius, and 
then saturates at k « 7.4 x 10^. The majority of pairs 
are produced at distances < R^s (see plot for N{< x)), 
which supports the assumption about the length of the 
cascade zone being ~ i?Ns made in § 3.1. The number 
of pair generations is highest at < i?NSi up to six in this 
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Figure 15. Cascade portrait of electron-positron pairs. Top panel, number of particles produced in each energy and distance bin dA^(a:, F) 
color coded in logarithmic scale according to the colorbar on the right. Middle left: E f{xi, E) - energy distribution of particles produced 
at four distances Xi^i = 1... 4; different line styles correspond to different distances according to the left plot legend. Bottom left: E f{E) 
- energy distribution of all particles produced in the cascade. Middle right: dN(x) — differential pair production rate - number of particles 
produced in a distance bin. Bottom right: N{< x) - total number of particles produced up to the distance x. Color lines in plots for 
E f{E), dN{x), and N{< x) show contributions of different cascade generations, lines are color-coded according to the right plot legend. 
Thick black lines show contributions of all cascade generations, x is the distance from the NS normalized to NS radius Hns and E is particle 
energy normalized to meC^. Particle number density is normalized to ncj. Parameters of this simulation: pulsar period P = 33 ms; the 
magnetic field in the PC has Bq = 10^^ G and pc = 8.8 x 10® cm; initial energy of primary particles = 2.3 x 10'. 


case®. The occurrence of most of the cascade generations 
at distances up to ^ i?Ns is the evidence that the photon 
mfp in a strong cascade is indeed small, ^ so that 
the cascade initiated by any given particle goes through 
several generations within the distance ~ i?Ns- The large 

® There are too few pairs produced in the 6th generation to show 
in the plot; curves corresponding to this generation are below the 
lower limit of all plots in Fig. 15; this generation shows in the 
portrait for photons, Fig. 16 


extent of the cascade zone relative to is due to con¬ 
tinuous injection of pair-producing CR photons (see the 
blue line in the plot of dN{x); generation 0 pairs are 
produced by CR photons). The largest contribution to 
pair multiplicity in this case comes from generation 1, 
pairs created by the first synchrotron photons (see plots 
for dN{x) and N{< x)). In our simulations for different 
pulsar parameters, contributions of generation 1 and 2 
to the pair multiplicity sometimes become comparable; 
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Figure 16. Cascade portrait of escaping photons for the same cascade as in Fig. 15. Notations and normalizations are the same as in 
Fig. 15. 


for weaker cascades the relative contribution of genera¬ 
tion 0 is higher than in this case, however we did not 
see generations 3 and higher producing the majority of 
pairs. The number of cascade generations is not very 
large in any of our simulations (several at most), but in 
each generations high numbers of pair producing photons 
are emitted which results in high multiplicity. 

The energy of created pairs decrease with distance (see 
plots for dN{x,E) and E f{xi,E)), mostly because of 
energy losses of primary particles which results in lower 
energy CR photons. The pair spectrum extends down 
to a few mc^, since the cascade is very efficient at con¬ 
verting initial pair energy into more photons (and pairs). 
Degradation of pair energies through cascade generations 
discussed in § 4.3 is clearly visible on the plot for E f{E) 


- the maximum pair energy systematically decreases with 
cascade generations at all distances. 

Figure 16 show the portrait of photons escaping the 
cascade. Photon generation 0 is CR while the higher 
generations (< 1) are synchrotron/cyclotron radiation. 
Although the highest energy photons are produced near¬ 
est the NS surface, these photons are absorbed by pair 
production attenuation so that the spectra at the lowest 
altitudes show sharp cutoffs near 100 mc^. This cutoff is 
clearly visible in the plot for Ef{xi,E) for xi = 0.15i?Ns 
shown by a dashed line. In the plot for dN{x,E) the 
cutoff is evident as a sharp horizontal boundary of the 
colored region for x < 0.16i?Ns- The escaping photon en¬ 
ergies increase with distance from the NS surface, as the 
magnetosphere becomes more transparent. The highest 


















































Pair multiplicity in young pulsars 


19 


escaping photon energies are produced near the end of 
the cascade, at around 4 NS radii. The highest energy 
photons escaping the cascade are CR photons (blue line 
in the plot for Ef{E)). Synchrotron radiation is emitted 
by pairs right after their creation, so that the pair forma¬ 
tion and synchrotron radiation end at the same distance 
- in the plot for dN{x), the number of photons in each 
generation drops at large distances in accordance with 
the drop of number of injected pairs shown on a similar 
plot in Fig. 15. Above ^ the emitted CR photon 
energies drop as the primary particles continue to lose 
energy. The spectrum of escaping photons also broadens 
at the lower end because, as the magnetic field decreases, 
so does the cyclotron energy which sets the lower limit of 
the synchrotron spectrum. Thus the lowest and the high¬ 
est energy escaping photons are produced at the largest 
distances from then NS, cf. spectra at different Xi in the 
plot for Ef{xi,E). The bulk of high energy emission 
from the PC cascade comes from synchrotron radiation 
of pairs in generations 1,2, and 3. 

Quantitatively, our semi-analytic theory compares 
with this particular numerical simulation as follows. For 
pulsar parameters used in this simulation, the energy 
of primary particles according to eq. (41) should be 
« 3.6 X 10^, which is « 1.6 times larger than the 
result obtained from numerical simulations with PAM- 
INA code « 2.3 x 10^. For the multiplicity of 

the CR-synchrotron cascade started by primary particles 
with monochromatic energies the semi-analytic 

model gives « 1.56 x 10®, from eq. (23), which 

is also about 2 times higher than the value obtained in 
numerical simulations « 7.4 x 10"*. The combined 
model from § 7, which uses the analytic model for parti¬ 
cle acceleration as an input for the semi-analytic model 
of CR-synchrotron cascade, predicts for the multiplicity 
« 2.9 X 10®, see Fig. 14, which is « 4 times larger than 
the multiplicity from numerical simulations. The dis¬ 
crepancy with the semi-analytic model for several other 
numerical simulations we performed with different pa¬ 
rameters is of the same order. We attribute this discrep¬ 
ancy mostly to the approximation of constant magnetic 
field - in the numerical simulations photon absorption de¬ 
creases and becomes less efficient with the distance. We 
think that for such a simple model the agreement with 
the numerical simulations is reasonable and the model 
can be used for estimates of multiplicities in young ener¬ 
getic pulsars. 

9. RESONANT INVERSE COMPTON SCATTERING 

Another emission mechanism for relativistic particles 
besides curvature and synchrotron radiation is inverse 
Compton scattering (ICS). In strong magnetic fields typ¬ 
ical for pulsar polar caps, ICS can occur in the resonant 
regime, when the photon energy in the electron’s rest 
frame is equal to the cyclotron energy. The cross-section 
for scattering of such photons is greatly enhanced com¬ 
pared to that of non-magnetic scattering. It has been 
noted that resonant ICS (RICS) with the soft thermal 
photons from the NS surface is important for high-energy 
emission from pulsar polar caps (Stumer 1995; Zhang & 
Harding 2000) and can be important in the development 
of polar cap cascades because for quite a wide range of 
pulsar parameters, scattered photons can be above pair- 


log AricS 



log e± 

Figure 17. Distance (in cm) over which a particle loses its energy 
via resonant inverse Compton scattering. Contours of logAj^j^g 
are plotted as a function of logarithms of particle energy e± and 
magnetic field B in Gauss for a NS surface temperature T = lO^K 
and fis = 0.5. 


formation threshold (e.g. Stumer et al. 1995; Zhang & 
Harding 2000). In this section we argue that although 
RICS can play a role in the development of polar cap cas¬ 
cades, it never becomes the dominant source for pair mul¬ 
tiplicity and, therefore, considering only CR-synchrotron 
cascades provides adequate estimates for pair multiplic¬ 
ity in strong cascades of normal pulsars. 

Eirst, let us consider the efficiency of RICS in trans¬ 
forming particle kinetic energy into radiation. The dis¬ 
tance over which a particle loses most of its energy to 
RICS is given by (Zhang & Harding 2000; Stumer 1995; 
Dermer 1990): 


Arics = -0.061 Tg 


In 


-1 


1 — exp — 


134Bi 


e± r6(l - fXs) 


cm (47) 


where Tg is the temperature of the NS surface in units 
of 10®K, Bi 2 is the magnetic field strength in units of 
10®^ G, and /is = cos 0s, where 9s is the angle between 
the momenta of the scattering photon and particle in the 
lab frame. If the NS is young, its surface temperature 
should be about 10®K, it emits X-ray photons necessary 
for RICS from the whole surface, and for particles at a 
distance comparable to the NS radius the range of /is is 
quite large. Eor cold NS, with full surface temperatures 
below fewx 10®K only the polar cap region, heated by the 
backflow of accelerated particles up to ^ 10®K, can emit 
enough photons for RICS to become important. In the 
latter case when the particle reaches a distance compa¬ 
rable to the width of the polar cap Vpc ~ 1.4 x 10^ P cm, 
where P is pulsar period in sec, the range of /is gets 
very small and photons quickly get out of resonance. So, 
for young NSs RICS is an important radiation process 
if Arics ~ ^ns; for old NSs this condition changes to 
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Figure 18. Diagram showing the general chain of physical pro¬ 
cesses in a strong polar cap cascade. Cascade generations are shown 
on the left - numbers connected by double arrows. In each genera¬ 
tion particles (electrons and/or positrons) produce photons (7 cr 

- via curvature radiation, 7syn - via synchrotron radiation, 7 r.ics 

- via Resonant inverse Compton scattering), which are turned into 
pairs of the next cascade generation. The CR-synchrotron cascades 
studied in detail in this paper are shown by solid arrow, dashed ar¬ 
rows show RICS initiated branches which are discussed only in § 9. 


AriCS ^pc 

In Fig. 17 we plot Arics as a function of particle en¬ 
ergy €± and magnetic field B,ioTT = 10®K and fig =0.5. 
It is evident from this plot that primary particles, with 
€± > 10®, lose a negligible amount of their energy via 
RICS, and curvature radiation is the dominant emission 
mechanism for generation 0 of strong cascades. For cold 
NSs, RICS of photons from heated polar caps might be 
an important emission mechanism only for a very nar¬ 
row energy range of low energetic particles, in one of 
the later cascade generation - the parameter space be¬ 
tween two Arics = 10^ contours is rather small - and for 
most values of B the scattered photons will be below pair 
formation threshold. Hence, RICS can be completely ne¬ 
glected in strong polar cap cascades of cold NSs. For hot 
NSs RICS becomes an important emission mechanism for 
a wide range of moderate particle energies e± < 10^ - the 
energy range between contours of Arics = 10® in Fig. 17 
is quite wide. In the latter case, the diagram for physical 
processes in a strong polar cap cascade can have the gen¬ 
eral form shown in Fig. 18, with RICS photons in some 
cases carrying non-negligible energy starting from cas¬ 
cade generation 1. The extension of the semi-analytical 
analysis of CR-synchrotron cascade developed in §§3-7 
to the whole cascade taking into account multiple cas¬ 
cade branches is not straightforward and we postpone 
it to future publications. Here we will make rough esti¬ 
mates of the contribution of RICS cascade branches to 
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Figure 19. Ratio of of the characteristic energy of RICS pho¬ 
tons RjQg emitted by the first generation pairs to the energy of 

escaping photons: contours of log are plotted as 

a function of of logarithms of the radius of curvature of magnetic 
field lines pc in cm and magnetic field strength B in Gauss. We 
used the following values for gap parameters P = 33 ms, = 2 
and Xa = 1/7. 



the multiplicity of the whole cascade. 

Let us now discuss the efficiency of RICS cascade 
branches in splitting the available energy into pairs. The 
energy of photons scattered by particles of energy e± in 
the RICS regime is (e.g. Zhang & Harding 2000) 

Cy.Rics = 2e±,F&, (48) 

where e±,F) given by eq. (12), is the kinetic energy of 
a particle moving along a magnetic field line - the final 
energy of freshly created pairs after they emit all per¬ 
pendicular to B energy via synchrotron radiation. Us¬ 
ing eqs. (48), (12) together with eq. (21) for the energy 
of CR photons and eq. (41) for the energy of the pri¬ 
mary particle, we get an upper limit on the energy of the 
generation 1 RICS photons - the highest energy RICS 
photons. In Fig. 19 we plot the ratio of the energy of 

generation 1 RICS photons to the energy of pho¬ 

tons escaping from the cascade zone Cesc as a function 
of the radius of curvature of magnetic field lines and the 
strength of the magnetic field B. It is easy to see that 
for magnetic fields weaker than a few xlO^^ G, even the 
highest energy RICS photons are not capable of produc¬ 
ing electron-positron pairs. For stronger magnetic fields, 
however, RICS photons do contribute to pair multiplic¬ 
ity. 

The characteristic energy of RICS photons in terms 
of the energy of the previous generation photon can be 
obtained by substituting e±_F into eq. (48) 

- 1/2 

(49) 
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log 

Figure 20. Ratio of the characteristic energies of synchrotron 
^ 7 ,'syn RIGS photons emitted by freshly cre¬ 

ated pairs ((i -(- 1)^^ generation cascade photons): contours of 
log |^^ 7 ,'syn/^^|,^~p\Qgj plotted as a function of logarithms of the 
(i) 

parent photon energy ^ and magnetic field strength B in Gauss 
for pc = lO^cm. 

This equation describes the energy degradation in each 
cascade generation for the RIGS process. In Fig. 20 we 
plot the ratio of the characteristic energies of synchrotron 
and RIGS photons (given by eqs. (49) and (16) corre¬ 
spondingly) produced by pairs created by the same par- 

(i) 

ent photon as a function of the parent photon energy ' 
and magnetic field strength B. The energy of RICS pho¬ 
tons are always smaller than the energy of synchrotron 
photons, and for B < 10^^ G significantly so. Because of 
the much faster energy degradation in the RICS process, 
cascade branches initiated by RICS photons are in gen¬ 
eral shorter than those initiated by synchrotron photons. 
In later cascade generations RICS photons freely escape 
the cascade zone while synchrotron photons emitted by 
the same particles are absorbed, still splitting the en¬ 
ergy of the parent photons into pairs. Therefore, RICS 
cascade branches should be in general less efficient in 
splitting the energy than synchrotron ones. 

Finally let us address the question of how much en¬ 
ergy is going into RICS branches of the cascade. The 
energy powering RICS branches is the kinetic energy 
of pairs moving along magnetic field lines, while syn¬ 
chrotron branches are powered by the perpendicular en¬ 
ergy of freshly created pairs. From eqs. (12), (13) we get 
for the ratio of the perpendicular to parallel pair energies 



Xa > b (see eq. (8)), the minimum value for this ratio 
is I — 0.414. In Fig. 21 we plot e_L/e|| as a function 
of the energy of the parent photon and magnetic field 





lOgCy 

Figure 21. Ratio of perpendicular to parallel to B energy of 
freshly created pairs: contours of ex/e|| as a function of logarithms 
of the parent photon energy and magnetic field strength B in 
Gauss for pc = lO^cm. 

strength B. For most of the parameter space in each 
pair creation event, the fraction of energy going into a the 
RICS cascade branch is smaller than those going into the 
synchrotron branch. Even in the case when more energy 
is left in the pairs’ parallel motion, the energy available to 
the RICS branch is only 1.5 times larger than the energy 
available for synchrotron branch. In this case, however, 
the value of Xa is close to 6, i.e. pair formation occurs 
near the kinematic threshold and is not very efficient (see 
§ 3.1). Hence, cascades where RICS branches have more 
available energy should not be very efficient. 

Summarizing the above arguments, RICS branches are 
less efficient in splitting photon energy into pairs and the 
energy available to those branches is in the best case com¬ 
parable to the energy available for synchrotron branches. 
The effect of RICS branches are strongest for higher val¬ 
ues of the magnetic field, > few x 10^^ G, when pair 
formation happens close to kinematic threshold and the 
cascade multiplicities are lower than for weaker magnetic 
fields. The final multiplicity of the total cascades should 
be less than twice that of the pure CR-synchrotron cas¬ 
cade in the best case and so the pure CR-synchrotron 
cascades studied in this paper provide good estimates 
for the multiplicity of strong polar cap cascades. 

10. FLUX OF PRIMARY PARTICLES AND PAIR 
YIELD 

As we have shown above, for the most efficient polar 
cap cascades the total pair yield for a given value of the 
magnetic field B depends mostly on the flux of primary 
particles. According to self-consistent models of polar 
cap acceleration zones (TIO, TA13), particle accelera¬ 
tion is intermittent and the pattern of plasma flow and 
acceleration efficiency depends on the ratio of the current 
density imposed by the magnetosphere to the GJ current 
density jm/jci as well as the boundary conditions at the 
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Figure 22. Snapshot of the phase space for Ruderman-Sutherland 
cascade at the end of a discharge. Cascade an pulsar parameters: 
Jm/jGJ = 1 B = 10^^ G, P = 33 ms, = 1. Horizontal axis 
- particle positron x, vertical axis — particle momentum normal¬ 
ized to rriec; the vertical axis is logarithmic except for the region 
around zero momentum (—5 < p < 5), where the scale is linear. 
The top panel shows phase space portrait of the cascade: each 
dot represents a numerical particle (every 10th particle is plotted); 
blue dots - electrons, red dots — positrons, black dots - photons. 
Three plots beneath show particle number density in phase space: 
p— vs X - electrons, vs a: - positrons, p^ vs x — photons. Par¬ 
ticle number density is color-coded according to the color map on 
the right in units of ncj. Particles which produce most of the pairs 
are positrons inside the area surrounded by dotted line marked as 
“bulk of pair-producing particles”. 

NS surface - whether particles can be extracted from the 
surface or not. 

There are essentially three qualitatively different 
regimes of plasma flow that determine the flux of primary 
pair-producing particles: (i) in space charge limited flow 
with 0 < jm/jcj < 1 particles are accelerated up to very 


Figure 23. Snapshot of the phase space for cascade in Space 
Charge Limited Flow regime. Cascade an pulsar parameters: 
jm/jGJ = 1.5 B = 10^^ G, P = 33 ms, = 0.25. Types of 
plots and notations are the same as in Fig. 22. 

low energies, e± ^ few, and no pairs are produced, (ii) in 
space charge limited flow with jm/jcj > 1 a significant 
flux of primary particles of the order of ~ Joj/e is accel¬ 
erated through the gap while the gap is moving toward 
the NS, (iii) in all other cases - i.e. for any current den¬ 
sity in Ruderman-Sutherland model and for jmjjai < 0 
in space charge limited flow regime - in each burst of pair 
formation a blob of primary particles is produced during 
formation of the gap, and then no significant amount of 
primary particles is created until the formation of the 
next gap. 

In Fig. 22 we plot as an example for case (iii) a phase 
portrait and densities of plasma and photons in the 
phase space for a Ruderman-Sutherland cascade with 
jm/jc] = 1- The process of gap formation for this case 
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is described at the beginning of §6.2. When the gap is 
formed, only a very few particles leak from the plasma 
tail - in this case electrons, visible as a line-like feature 
in the top two panels of Fig. 22 - are accelerated in the 
gap, the flux of these particles is much less than ^ jcj/e 
and their contribution to the pair production can be ne¬ 
glected. The blob of primary particles which was created 
during the gap formation has a density of a few Uqj in its 
densest parts and a size comparable to the gap’s height. 
Those particles - in this case positrons in the area sur¬ 
rounded by a dotted line on the plot for positron phase 
space density in Fig. 22 - create the vast majority of 
pairs. The density of primary particles in the blob is 
high but no new primary particles are created until the 
next gap is formed, the duty cycle of the cascade is small. 
The resulting flux of primary particles averaged over the 
period of gap formation is rather low. 

In case (ii), in space-charge limited flow with super GJ 
current density, the situation is qualitatively different in 
that during the lifetime of the gap, and not only during 
the time the gap is forming, as in case (hi), a significant 
constant flux of particles is going through the accelera¬ 
tion zone. The duty cycle of such a cascade should be 
substantially higher than in any other types of cascades. 
As an example of such a cascade we plot in Fig. 23 a 
phase portrait and the densities of plasma and photons 
in phase space for a cascade in space-charge limited flow 
regime with jm/jci = 1.5. The gap is formed at some dis¬ 
tance from the NS and moves toward it. When the first 
particles are formed the process of electric field screen¬ 
ing proceeds in a similar way to the case of Ruderman- 
Sutherland cascades; the blob of ultrarelativistic parti¬ 
cles is created and particles leaking from it create a tail 
of mildly relativistic plasma screening the electric field 
behind the blob; this blog moves toward the NS creat¬ 
ing pairs. However a constant flux of particles extracted 
from the NS surface - in this case electrons, visible as 
a line-like feature in the top two panels of Fig. 23 - are 
accelerated in the gap. The density of these electrons 
is high (for the case shown in Fig. 23 it is ~ 1.25nGj) 
and these particles, and not the particles from the blob, 
produce most of the pairs. 

Because of the intermittency of pair formation, the 
resulting multiplicity must be adjusted by the relative 
fraction of time during which primary particles are pro¬ 
duced. If Tactive is the time of active particle acceleration 
and Tcascade is the time between the beginning of suc¬ 
cessive bursts of pair creation, then the pair multiplicity 
as discussed in previous sections should be multiplied by 
the attenuation factor 

n "^active /r"i\ 

Ik = Tf, - (51) 

to get the average multiplicity of pair cascades. 

The existing self-consistent simulations of particle ac¬ 
celeration (TIO, TA13) are inconclusive about the cas¬ 
cade repetition rate. The numerical resolution in these 
simulations was inadequate for that purpose - because 
of copious pair formation the Debye length of plasma 
at some point became smaller than the cell size and the 
formation of the plasma tail could not be simulated accu¬ 
rately; the repetition rate, however, is determined by the 
mildly relativistic plasma in the plasma tail, as the next 
burst of pair formation starts only when plasma leaves 


the polar cap region. Future numerical simulations would 
address this issue and provide accurate values for cascade 
duty cycles; meanwhile here we will try to give a rough 
estimate for the cascade repetition rate using the follow¬ 
ing simple physical arguments. 

The strongest pair formation can continue up to the 
distance of the order of a few NS radii i?Ns- Although 
the exact physical mechanism responsible for formation 
of the mildly relativistic plasma tail is not clear (see pre¬ 
vious paragraph) it seems to us plausible to expect that 
when plasma injection due to pair formation stops, rever¬ 
sal of some pairs toward the NS - plasma leakage from 
the blob to the tail - should taper off as well. So, the 
time when plasma is filling the polar cap region is of the 
order of ~ i?Ns/c- The mildly relativistic pairs which 
fill the polar cap region after the blob of plasma has 
moved away are moving with relativistic velocities, so 
the time the plasma needs to clear out is also of the or¬ 
der of ^ Rns/c. Therefore, the cascade repetition rate 
seems to be Tcascade ~ few x Rns/c 

For the case (iii), in the flow regime with vacuum gaps 
- Ruderman-Sutherland cascades or space-charge limited 
flow in return current regions - particles are accelerated 
only during formation of the gap, the time of active par¬ 
ticle acceleration is quite short. The longitudinal size of 
the blob of primary particles is about the size of the gap 
/igap, and the primary particle density in the blob is a few 
ricj. All primary particles in the blob move in the same 
direction and the blob passes any given surface normal 
to the field lines during time Tactive ~ ^gap/c- If the next 
gap will form after the blob of primary particles moves a 
few Rfis away, the attenuation factor should be 

/. ~ (p = 4 X 10-= (52) 

^NS 

For the Crab pulsar this attenuation factor is around 
0.001 if we assume P = 33 ms, P = 3 x 10^^ G, Xa = 
1/7, and = I (for inclination angle 60° and current 
density jm — jcj)- The resulting multiplicity of such 
cascades in Crab would be < IO^Ugj- 

In the space-charge limited flow regime with super-GJ 
current density, case (ii), the cascade efhciency should 
be much higher. The gap appears at some distance and 
moves toward the NS (see § 6.2 in TA13). The process 
of gap formation is very similar to that in cascades for 
case (iii) discussed above. A blob of relativistic particles 
is created, which in this case propagates toward the NS, 
and pair plasma leaking from it creates a plasma tail fill¬ 
ing the region behind the blob which prevents formation 
of the next gap. But in this case particles extracted from 
the NS surface flow into the magnetosphere through the 
gap and are accelerated all the time the gap exists; these 
particles create most of the pairs. When the gap disap¬ 
pears, particle acceleration stops. The remaining plasma 
should leave the gap and flow into the magnetosphere be¬ 
fore the next gap can appear. The time necessary for the 
remaining plasma to leave should be comparable to the 
time of the gap’s existence, as all particles are relativis¬ 
tic. Therefore, in this case the attenuation factor should 
be /k I/few - much higher than given by eq. (52). So, 
polar cap regions with super-GJ current density should 
have a pair yield of about lO^ncj- 

Let us summarize our findings. Pulsars with polar cap 
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cascades operating in the Ruderman-Sutherland regime 
would not be efficient pair producer, the highest yield 
should be less that ^ lO^ncj- If the polar cap cascades 
operate in the space charge limited flow regime, then in 
regions with super-GJ current densities the pair yield 
is quite high, around lO^noj; in the regions with return 
current (anti-GJ current density) the yield should be less 
that ~ lO^riGj) and in regions with sub-GJ current densi¬ 
ties no pair plasma would be produced at all. We there¬ 
fore conclude that the maximum pair yield of a young 
pulsar would be less that lO^ncj, as only a fraction of 
the polar cap can have super-GJ current density. 

Pair yield is mostly affected by the duty cycle of the 
cascade, and not by pulsar parameters such as magnetic 
field strength B, radius of curvature of magnetic field 
fines pct or pulsar period P. This holds only for the 
regime in which the gap height is much smaller than the 
PG radius, true for all young pulsars. When the gap 
height approaches and exceeds the PC radius, the pair 
yield drops quickly, as shown in Figure 9. The distri¬ 
bution of jm/jci in the polar cap is determined by the 
pulsar inclination angle and we expect that inclination 
angle should be the most important factor determining 
pair yield of a young pulsar. 

11. DISCUSSION 

We have performed a systematic study of electron- 
positron pair cascades above pulsar polar caps for a vari¬ 
ety of input parameters including surface magnetic field, 
pulsar rotation period, primary particle energy and mag¬ 
netic field radius of curvature. We have also studied here 
for the first time the multiplicity of self-consistent pair 
cascades, i.e. those that are capable of generating cur¬ 
rents consistent with global magnetosphere models. 

We find that pair multiplicity is maximized for a mag¬ 
netic field strength near 10^^ G, independent of the other 
parameters. This value of field strength strikes a balance 
between maximizing the fraction of photons that pair 
produce - the weaker the magnetic field the higher the 
number of escaping photons - and the fraction of pho¬ 
ton energy going into pair production - in stronger fields 
photons are absorbed when they have smaller angles to 
B, and created pairs, having smaller perpendicular to B 
momenta, emit less synchrotron photons, which produce 
the next generation of pairs. This is true for curvature 
radiation - synchrotron cascades, which should be the 
dominant source of pairs in young energetic pulsars. For 
stronger magnetic fields B >3 x 10^^ G resonant inverse 
Gompton scattering of soft X-ray photons emitted by the 
NS surface can tap some of the pair’s energy parallel to 
B and increase the cascade multiplicity. RIGS cascade 
branches, however, are less efficient than the synchrotron 
branches; this together with the decrease of the photon 
absorption cross-section for near threshold pair creation 
in strong magnetic fields should not change the fact the 
polar cap cascade multiplicity reaches its maximum for 
B ~ 1012 G. 

We find that the pair multiplicity at the peak of the 
cascade cycle, k ^ 10®, is remarkably insensitive to pul¬ 
sar period, magnetic field and radius of curvature of mag¬ 
netic field lines. The reason for this is self-regulation of 
the accelerator by pair creation: for pulsar parameters 
resulting in more efficient pair production, the size of 
the acceleration gap is smaller and the primary particle 


energy is lower and vice-versa. The most important fac¬ 
tor determining the multiplicity of polar cap cascades is 
the flux of primary particles which depends on the “duty 
cycle” of the particle acceleration in time-dependent cas¬ 
cades. Estimating the “duty cycle” we find that the time- 
averaged pair multiplicity is limited to ^ 10® for the 
case of space charge-limited flow (free particle extraction 
from the NS) with super Goldreich-Julian current density 
(j/J gj > 1) and to only ~ 10® for the case of Ruderman- 
Sutherland gaps (no particle extraction from the PG) 
and space charge-limited flow with anti Goldreich-Julian 
current density {j/jcj < 0). The current density dis¬ 
tribution in the PG depends on the pulsar inclination 
angle, and so for young energetic pulsars the inclination 
angle should be the most important factor determining 
the total multiplicity of plasma in the magnetosphere. 

Our finding of a peak multiplicity from self-consistent 
pair cascades in young pulsars (^ 10®), that is largely 
independent of their period and magnetic fields, is very 
different from previous predictions from steady cascades 
that produce a range of multiplicities (10^ — 10'^) that 
depend strongly on pulsar parameters. This difference 
primarily results from much smaller gap heights in the 
time-dependent cascades, that keep the gap size much 
less than the PG radius over a large parameter range. 
This in turn results in a nearly uniform magnetic field 
strength throughout the gaps, higher accelerating elec¬ 
tric fields, higher primary particle energies and more ef¬ 
ficient pair production - primary particles are accelerated 
faster and photons are injected in the region of stronger 
magnetic field than in previous steady cascade models. 

We have also numerically simulated the photon spec¬ 
tra from self-consistent polar cap pair cascades. Our re¬ 
sults show that the cascade photon spectra from the more 
compact gaps nearer to the neutron star surface have cut¬ 
offs that are around 10 - 100 MeV. This is at the lower 
end the Fermi energy band (30 MeV - 300 GeV), which 
may explain why Fermi has not seen strong evidence for 
7 -ray emission from the PG. High energy emission from 
pair cascades is thus expected to occur in lower energy 
bands, below 100 MeV. PC cascade emission may have 
been detected recently from PSR J1813-1246 (Marelli 
et al. 2014) at X-ray energies by XMM-Newton and 
Chandra. The unusual X-ray light curve that shows two 
peaks separated by 0.5 in phase but offset by a quar¬ 
ter of a period with the y-ray peaks can be explained 
through a geometric model placing the y-ray emission in 
the outer magnetosphere and the X-ray emission at lower 
altitude above the PCs, and could be coming from the 
PC cascades. 

Our results show that PC cascades are more efficient 
in producing electron-positron pairs that was previously 
assumed. However, even the higher pair multiplicity of 
10® is not enough to account for the fluxes of parti¬ 
cles needed to explain the synchrotron radiation from 
PWNe. Based on a multiplicity of 10®, we estimate that 
the Crab pulsar produces a pair flux from each PC of 
about 2 X 10®® pairs s“®. The flux from both PCs is an 
order of magnitude smaller than the pair flux required to 
account for the radiation from the nebula, which is esti¬ 
mated to be about ^ 4 x 10"*® pairs s“® (de Jager et al. 
1996). Cascades in the outer magnetosphere are not very 
efficient pair producers (e.g. Hirotani 2006), and so the 
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injection of plasma by pulsars can not account for the 
population of particles in PWNe emitting at radio wave¬ 
lengths. These radio emitting particles must then have a 
different origin from particles emitting at shorter wave¬ 
length; for example, they might be picked up from the 
gas filaments in the supernova remnant or be remnants 
of some unknown acceleration mechanism in the early 
history of the nebula (Atoyan & Aharonian 1996). 

Another implication of our results concerns observa¬ 
tions of cosmic-ray electrons and positrons at Earth that 
have shown an excess of positrons over what can be 
produced in secondary cosmic-ray interactions (Adriani 
et al. 2009; Accardo et al. 2014), indicating the exis¬ 


tence of primary positron sources in the Galaxy. Var¬ 
ious studies (e.g. Gendelev et al. 2010) have estimated 
that PWNe could account for the excess of cosmic ray 
positrons. With higher pair multiplicity from young pul¬ 
sars, and subsequent acceleration of the pairs at the pul¬ 
sar wind termination shock, PWNe could produce a more 
significant contribution of primary cosmic ray positrons. 

This work was supported by a NASA Astrophysics 
Theory grant and a Fermi Cycle-5 guest investigator 
grant. 


APPENDIX 

A. ALGORITHMS FOR SEMI-ANALYTICAL CALCULATION OF CASCADE MULTIPLICITY 

Here we show pseudo-codes of algorithms used to compute cascade multiplicity. For calculation of Xai^j, B, pc) we 
computed and stored a table of 1 /xa values for a uniformly divided grid 77 x 30 x 20 in log x log B x log Pc space, 
and then used cubic piece-polynomial interpolation to get Xa for parameter values required by expressions used in the 
algorithms. 

Algorithm 1 computes total number of particles produced in synchrotron cascade initiated by a single primary photon 
with the energy Cj, physical processes are described in §4.3. e_phot and n_phot are the number and the energy of 
synchrotron photons emitted by particles of the current generation, e_phot is the energy of escaping photons, and n 
is the total number of particles produced in all previous generations. 


Data: e-y - energy of primary photon, Sesc - mfp of escaping photon 

Result: A^syn - total number of particles produced in synchrotron cascade initiated by a photon with the energy e-y 


Function NsynU-/, Sesc)'- 

C-CSC <— €7,esc(sesc) 
e_phot <— e-y 
n_phot <— 1 
n <— 0 

while e_phot > e_esc do 
n 4— n + 2 n_phot 
n_phot <— n_phot nsyn(e-phot) 
e_phot <— e 7 ^^\e_phot) 

end 

return n 


// eq. (10) 


// eq. (15) 
// eq. (16) 


Algorithm 1: Multiplicity of synchrotron cascade 


Algorithm 2 computes the total multiplicity k of a CR-synchrotron cascade according to eq. (23). e_part is the 
energy of the primary particle at the current distance sp, e_phot is the energy and n_phot is the number of CR photons 
emitted by the primary particle at the current distance, Nsyn is the number of particles generated in synchrotron 
cascades initiated by the primary particle at the current distance. Integration over the distance is done with the 
simple trapezoidal rule. 


B. TRANSITION TO RADIATION REACTION LIMITED REGIME 


Radiation reaction begins playing an important role in the dynamics of particle acceleration if particle energy losses 
per unit time become comparable to the work done on the particle by the accelerating electric field. For curvature 
radiation the condition for applicability of the free acceleration regime is (see eq. (17)) 


2 



(HI) 


The electric field E grows linearly with the distance l± traveled by the particle in the gap, eq. (27); the particle energy 
€± increases as eq. (32), and at some distance the condition (BI) is violated. However, if condition (Bl) holds 
up to the distance Zj. where the particle emits “gap-terminating” photons, then the gap length and the final particle 
energies can be well described by the free acceleration regime, neglecting radiation energy losses. Using eq. (27) for E, 
eq. (32) for e±, and substituting Z±,gap fo'' from the eq. (39) inequality (BI) becomes the condition on the magnetic 
field strength 

B > \afXaBq ~ 2.8 x IQ-^R, = 1.15 x lOi^G. 
o 


(B2) 
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Data: - energy of the primary particle, scr ~ characteristic size of the cascade, Sesc - nifp of escaping photons 

Result: total number of pairs produced by the primary particle with the energy 


begin 

divide [0, scr] in subintervals sq, •.., 

Nsyn.last 0 

n 0 

for ^ 1 to imax do 

// CR radiation at distance Sj 
e.part ^ e±(e![, Si) 
e_phot e-y, cR(e-part) 
n_phot ■<— ncR(e_part) 

// synchrotron cascade multiplicity 
Nsyn <— n.phot A^syn (e-phot, Sesc) 

// trapezoidal rule 
n -4— n + 0.5(Nsyn + NsynJast)(sj — Si_i) 
Nsyn.last ■<— Nsyn 

end 

return n 

end 


// eq. (19) 
// eq. (21) 
// eq. (22) 

// algorithm 1 


Algorithm 2: Multiplicity of CR-synchrotron cascade 


In the second step we used the value Xa = 1/7- The radius of curvature pc and pulsar period P cancel out, except for 
a weak dependence of Xa on pc- We see that for most non-millisecond pulsars radiation reaction can be neglected. 


C. LIMIT ON 1-D APPROXIMATION 


The one-dimensional approximation for particle acceleration in the gaps works well if the length of the gap is much 
smaller that the width of the polar cap. As the length of the gap is ^ 2 times larger than Zgap, we can set the formal 
limit on 1-D approximation as 

^gap ^ I"pc 5 (^7) 


where rpc = Pc is the polar cap radius. Using eq. (37), (39) for /gap and we get the following limits on 

pulsar parameters for applicability of the 1-D approximation 


^4/7p-13/14 > 7 / C y/" / 

6 \2ttRIJ ^ 


1/7 




(C2) 


The pulsar magnetic field is usually estimated assuming magnetodipolar energy losses from the values of period P and 
period derivative P as B = 3.2 x 10^® \/ PP. Expressing B in this way and putting the numerical values for physical 
constants, we get from eq. (C2) the following condition for applicability of 1-D approximation in terms of P and P 

P > 3.21 X 10-22;yy2^-3/2^^ p9/4 _ 4 3 ^ ^ p9/4 ^ 

in the last step we used values = 2 and Xa = 1/7. For cascades along dipole magnetic field lines at the edge of the 
polar cap with the radius of curvature 

Pc = ~ 9.2 X lO'^yPcm, (C4) 

^ ^pc 


where 0 pc 


y^27ri?Ns/cP is the colatitude of the polar cap edge, condition (C3) takes the form 


P > 4 X lO-i^P"/"^. 


(C5) 
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